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SYNOPSIS 


P.K. RAIMA 
Ph.D, (Physics) 

Indian Institute of Technology Kanpur 
FEBRUARY 1987 

FOm FACTORS AND TRANSITION CHARGE DENSITIES FOR 
THE MULTIPCLE ELECTROEXCITATION OF SOME NUCXEI 
IN THE MASS 40 TO 70 RB3ION 

Inelastic electron scattering has proved to be an 
excellent method for exploring nuclear structure. A measure- 
ment of the electroexcitation cross section permits one to 
obtain nuclear dynamical properties such as the transition 
charge and current densities. The calculation of these 
transition densities in terms of the nuclear wav ef unctions 
provides a sensitive test of the latter. 

Recent inelastic electron scattering experiments 
involving a number of doubly even 2p-1f shell nuclei have 
provided valuable data for the quadrupole and hexadecupole___ 
electroexcitation of the yrast levels with J =2 and 4 in 
these nuclei. Within a measured momentum- transfer range upto 
3 fm“^ the quadrupole form factors display two distinct 
maxima appearing at q 0,7 and 1,6 fm“^ , The hexadecupole 
form factors, on the other hand, are characterized by maxima 
occurring at q 1 ,2 and 2,2 fm~^ , The relative magnitude irf 
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the form factors at the two maxima is expected to provide 
a sensitive test of the nuclear microscopic models involved 
in the calculations. 

In recent years considerable attention has also been 
devoted to the coordinate— space reconstruction of the trans.i— 
tion charge densities associated with the electroexcitation 
of nuclear levels; this has been achieved by employing the 
Fourier Bessel transforms in conjunction with the distorte(^. 
wave Bom-approximation, The empirical transition charge den- 
sities acquire additional significance vis— ci-vis the test of 
the model wave functions since they contain structural infor- 
mation unaffected by the accuracy of the usual plane-wave 
Born-approximation for large momentum transfers. 

In a recent shell-model study the quadrupple electro- 
excitation form factors as well as. the transition charge 
densities for some doubly even Ni isotopes were calculated in 
terms of the shell model wave functions resulting from semi- 
empirical effective interactions operating in the restricted 
valence spaces involving the [(2p2y2» ^^5/2* 

(If^/g)*"^ (2P3/2. configurations. The use 

of the shell model wave functions involving limited number ef 
configurations, however, resulted in a highly mass-dependent 
set of proton and neutron effective charges. 

The available form factor data in doubly even Zn 
isotopes (with A 66) has not been analysed in a microscopic 



AJLJ.X J 


framework thus far, A conventional shell model description 
of these isotopes in terms of just one 

expected to be quite inadequate if one considers the fact 
that the available data on the reduced E2 transition probabi- 
lities suggests enhanced rotational collectivity compared to 
that in the Ni isotopes. In view of the intractability of 
a realistic shell model calculations it is desirable to have 
an approximate yet reliable ansatz for the multiparticle wave 
functions generated within the full 2p-1f valence space; 

In the present work v/e show that the yrast wave functions 

It 4 . 

for the J = 0 , 2 levels projected from the Hartree-Fpck— 

Bogoliubov (HFB) intrinsic states resulting from realistic 
effective interactions operating in the full 2|>-tf shell 
permit a fairly satisfactory description of the 2"*" -form factors 
in some doubly even Ni and 2n isotopes (with 58 < A _< 68), 

In this context wo also, discuss the calculation of the transi- 
tion charge densities in terms of the self-consistent HFB states. 
The calculations dem.onstrate that the use of the projected HFB 
wave functions obviate the necessity of employing widely diff- 
erent effective charges for various nuclei. 

In Chapter I we briefly outline the overall perspective 
of the calculations presented in the thesis. 

In Chapter II we discuss in some detail the use of the 
projected HFB v^ave functions vis-^~vis the calculation of the 
form factors as well as transition charge densities associated 
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with the electroexcitation of yxast levels. 

We next discuss in Chapter III the results for the 

form factors, transition densities as \\fell as the E2 transi- 

+ + 

tion strengths for the 0 -* 2 transition in some doubly 
even Ni and Zn isotopes. It is seen that a comparison of 
the projected HFB transition densities with the ones extra- 
cted from the observed form factor data (via the Fourier- 

Bessel analyses) offers important guidelines vis-l-vis the 

40 

choice of a model of the Ca core-polarization contri- 
butions , 

In Chapter IV we have tested the efficacy of the 
projected HFB description vis-S-vis the recently-measured 

Tt + 

data for the electroexcitation of J =4 levels in a 
number of dcubly-even isotopes of Ti,Cr,Fe and Mi, Some 
recent shell model calculations involving restricted configu- 
ration mixing have required highly mass-dependent effective 
charges (for the 4 electroexcitation) viiich are also drasti- 
cally different from the set employed for a reasonable 
quantitative agreement with the observed ( 0 ■* 2 ) form 
factors in these nuclei. In contrast to this, our calculations 
provide a satisfactory quantitative discussion of the 
( 0 4 ) form factors in terms of the effective charges quite 

close to the ones required for a ccmparable fit to the 
(0^-» 2 "*") data in these nuclei. 

The recent measurements of the form factors involving 
+ + 

the 0 4 transitions has also provided reasonably accurate 



estimates of the heretofore unmeasured E4 transition proba- 
bilities, B ( E4f o"^ •* 4 "^) via an extension of the data_ into 
the momentumr-transf er of the photon point. In this context 
we demonstrate in Chai^ter IV the quantitative as well as 
qualitative unreliability of the Bohr and Mottelson 
prescriptions for relating the E4 transition matrix elements 
with the hexadecupole intrinsic moments. 

Finally in Chapter V we summarize the results and 
discuss some possible extensions of the present work. 



CHAPTH^ I 


IMTRODUCTION 


The elastic as well as inelastic el ectron~nucleus 

scattering turns out to be an excellent tool for studying 

1 2 

nuclear structure for a number of reasons 

One of the important features that characterizes 
the choice of the electron as a probe is the weakness 
of the electromagnetic force responsible for the electron- 
nucleus interaction. Since the coupling constant for the 
electromagnetic interaction is much smaller than the 
strength of the nuclear interaction, the electron hardly 
perturbs the nucleus. The usefulness of the electron 
as a probe is also augmented by the fact that the inter- 
action between the electron and the nucleus is the well- 
understood electromagnetic interaction that is described 
by an exact theory—" quantum electrodynamics. The 
nuclear properties can therefore be extracted from the 
observed data in an unambiguous^ quantitative manner. 

An important consequence of the weakness of the 
electromagnetic interaction is the possibility of providing 
an adequate quantitative discussion of the electron- nucleus 
scattering in terms of the first Born approximation (invol- 
ving only the one-photon exchange graphs) in conjunction 
v/ith the description of both the incoming and outgoing 
electrons as plane waves. The single-step nature of the 
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scattering process underlying this approximation implies 
that a given momentum transfer is a direct measure of the 
spatial resolution of the probe. 

In the context of the excitation of nuclear levels* 

the inelastic electron scattering has advantages over the 

processes involving real photons ; for a fixed energy loss 

of the electron one can vary the momentum q transferred 

to the nucleus. This contrasts keenly with the situation 

in the case of the real photons where one has only a single 

possible momentum, transfer for a given energy transfer. Thus 

inelastic electron scattering facilitates the study of the 
2 

complete q behaviour of the matrix elements involved in 
the transition. 

In the framework of some phenomenological models for 
inelastic electron scattering, the differential cross-section 
can be characterized approximately by , v\here Is 

the spherical Bessel function, L is the multipolarity of 
the transition and R the nuclear dimensions. As one increases 
q, by either raising the incident energy or by considering 
large-angle scattering, the factor qR) approaches its 
first diffraction minimum located, say, around q^ , However, 
higher multipoles (L+l) are expected to rise up to their 
first maximum for q ft*q^.. One may thus maximize a certain 
multipole while suppressing the strength of lower, multi poles. 
This raises the possibility of exciting, strongly new levels 
which would not be excited by other means; 



We must also mention here some of the inherent draw- 
backs associated with the electro-excitation experiments^ 

The significant distortion of the electron wavef unction 
by the Coulomb interaction necessitates the use of distorted 
wave functions for interpreting experiments involving hesfvy 
nuclei (2 > 20), Further, due to their small mass, the 
electrons are easily deflected. The resulting radiation 
provides strong background to the scattering cross sections^ 

In the case of nuclear electro-excitations, one can 
extract the Fourier transforms of the charge and current 
densities associated with the transition, from a measurement 
of the cross section by invoking the plane wave Born appro* 
ximation. The momentum- transfer dependence of the nuclear 
matrix elements contains information about the spatial 
structure of the nuclear states, A comparison of the observed 

momentum-transfer dependence with the one resulting from 

/ 

model wavefuncti ons is expected to provid^'sigriif leant test 
of the latter, 

A number of recent inelastic electron scattering 

*1* Hh 

experiments have provided valuable data on the 0 •* 2^ 2 

as well as 4| transitions in several 2p-1f shell 

3-1 2 

nuclei , The available form factor data covering the 
range q = 0,1 - 3,0 fm*"^ is characterized by the following 
salient features: 

( i) The observed form factors for the transition ^ 2.^ 2 

display two distinct maxima appearing at q 0,7 and 1,6 fm ~^4 



The minima of the form factors appear at q ~ 1 ,3 fm 

( ii) The maxima of the form factors associated with 
the o’*" *• 4^ transition occur at q '^1,2 and 2,2 fm”^ , 

( iii) Although the qualitative features of the form 
factors associated with the 0^ •* 2^ transition are quite 
similar to those of the observed form factors for the 

4- 4. 

0 2 ^ transition, the magnitudes in the case of former 

transitions are smaller by an order of magnitude, 

(iv) In the case of the ** 2| transitions, the magnitude 
of the form factors at the second maxima are usually 10 to 
15 times smaller compared to the values at the first maxima. 
In the case of the 0 •• 4^ transition, on the other hand, 
the ratio of the form factors at the first and the second 
maxima ranges between 50 and 100; The ratio of the form 
factors at the first, and second maxima is expected to , 
depend sensitively on the nuclear microscopic model Involved 
in the calculations; 

The topic of the coordinate space reconstruction of 

the transition charge densities associated with the electrort 

excitation of nuclear levels has received a lot of attention 
1 2 

in recent years * The transition charge densities have 
been deduced from the available inelastic scattering data 
in a nearly model-independent way by employing the technique 
of the Fourier-Bessel transforms in conjunction with the 
distorted-wave Born approximation. If one considers a simple 


single- particle transition ““ an excitation of a nucleon 
from state a to state b — the transition charge density 
is essentially the product of the radial v/avefunctions of 
the initial orbit a and the final orbit b times the charge 
of the nucleon. The empirical transition charge densities 
thus provide a detailed picture of the underlying nuclear 
wavefunctions , Further, they acquire additional signifi- 
cance vis— a— vis the test of the model wavefunctions since 
they provide structural details unaffected by the accuracy 
of the usual plane— wave Born approxiriHtion for large momentum 
transfers. 

In this thesis we attempt a consistent microscopic 
description of the available experimental data on the form 
factors as well as transition charge densities for the 
2p-1f shell nuclei in terms of the angular momentum projected 
Hartree-Fock-Bogoliubov (HFB) wave functions^ ^ resulting 
from the realistic effective interaction^^ obtained by Kuo 
and Brown for the valence space spanned by the ^^ 3 / 2 * 

^^1/2* ^^5/2^ orbits. In view of the intractability of the 
unrestricted shell model calculations involving the 

^p3/2» ^^1/2* ^^5/2^ * space for the nuclei with 40<A<80 
and this is because the dimensionalities of the relevant 
shell model matrices are prohibitively large — the approach 
discussed in the present work is expected to provide an 
approximate yet reliable ansatz for the multiparticle wave 
functions generated within the full 2p-1f valence space,. 



In Chapter II v/e describe the calculational framework 
involving the use of the projected HFB wave functions for 
interpreting the recent data on the form factors as well 
as transition charge densities associated with the electro- 

+ 4 - 

excitation of the 2^ and 4^ levels. 

We next discuss in Chapter III the results for the 
form factors, transition charge densities as well as the 
E2 transition probabilities for the 0*^ -*> 2^ transition in 
a number of doubly even Ni and Zn isotopes in the second 
half of the 2p-1f shelly 

In a recent study ‘ , the shell model wave functions 
resulting from semi- empirical effective interactions operat- 
ing in the restricted valence spaces involving the 

l-^^P3/2» ^^5/2’ ^^1/2^ ^^^7/2^ ^^^^3/2^ ^^5/2* ^^1/2^'^^”* ^ 

-f 4 

configurations have been employed for calculating the 0 •* 2^ 

form factors as well as the relevant transition charge 
densities for some doubly even Mi isotopes. The inadequacy 
of these shell model wave functions for interpreting and 
correlating the available data is reflected in the use of 
significantly mass-dependent proton and neutron effective 
charges that are required for optimizing the agreement with 
the experiments. 

In the case of the doubly even Zn isotopes, the 
inappropriateness of the shell model descriptions involving 

It V 

just one (Ify/o) * hole is emphasized by the fact that the 
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available data on the reduced E2 transition probabilities 

indicates considerably enhanced rotational collectivity 

18 19 

compared to that in the Ni isotopes * , Because of these 

considerations, the available form factor data in doubly 
even Zn isotopes (with A 66) has not been interpreted 
in a microscopic framework thus far. 

The calculations presented in Chapter III demons- 

TC 

trate that the yiast wave functions with J = 0 , 2 

projected from the HFB intrinsic states resulting from 

realistic two-body eff active interactions operating in the 

full 2p-1f shell yield a reasonably satisfactory unified 

description of the available form factor data for the 

o'*' 2^ transition in the nuclei and ^ 

It turns out that the use of the projected HFB wave functions 

obviates the necessity of employing videly different effective 

charges for various nuclei, V/e have also discussed the 

calculation of the transition charge densities in terms of 

the self-consistent projected HFB wave functions. It is 

seen that a comparison of the model space transition densities 

with the empirical ones extracted from the data via the 

Fourier-Bessel analyses supports unambiguously one of the 

40 

two existing models for the Ca core-polarization transition 
density. 

The calculations presented in Chapter IV examine the 
efficacy of the projected HFB v^ave functions vis-^-vis the 
recently-measured data involving the hexadecupelo 



electroexci tation in a number of doubly even isotopes of 
Ti, Cr, Fe and Ni, It is encouraging to note that, in 
the framework cf the projected HFB method, the effective 
charges required for a satisfactory quantitative discussion 
of the 0 ■* 4^ form factors are c^ite close to the ones 

required for a comparable fit to the quadrupole electro- 
excitation data in these nuclei. This again contrasts 
keenly with the results obtained in some recent shell model 

17 

calculations involving restricted configuration mixing 

in which the optimization of the agreement vdth the 

observed IF 1 (0 -♦ 4 ) data has required highly mass - 

dependent effective charges which are also drastically 

different from the values employed in the (shell model) 

“i" 

studies of the 0 -* 2^ form factors. 

An important implication of the recent measurements 

+ + 

of the form factors associated with the 0 4, transitions 

1 

is the availability of the heretofore unmeasured E4 transition 
probabilities, B( E4; O' 4 ), via an extension of the data 
into the momentum-transfer of the photon point. In view 
of this, we digress a little in Chapter IV to demonstrate 
the quantitative as well as qualitative unreliability of 
the usual Bohr and Mottelson prescriptions for relating the 
E4 HBtrix elements and the intrinsic hexadecupole moments. 

It turns out that the relative efficacy of the rotor model 
predictions for the quadrupole and hexadecupole operators 
is mostly dictated by the degree of closeness between the 



eigenfunctions of these operators and the ones resulting 
from the two-body interactions via the HF/HFB prescrip- 
tion. 

Finally in Chapter V we summarize the main results 
and discuss some possible extensions of the present work. 
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CHAPTER II 


MICROSCOPIC DESCRIPTION CF THE . ELECTROE^XITATION - 

'T!rTmrs“”(5FT«rp'ffs?e^£c ^I'Kge-PdciK-BciG’aTossv"'' 

Method ^ ^ 


11.1 Introduction 

In this chapter we shall first present in Section 

11.2 the theoretical framework relevant for a general 
discussion of the nuclear electroexcitation in the context 
of the plane~wave-Bori>-approximation (PWBA), The informa- 
tion about the nuclear vrave function is contained in the 
(nuclear) form factors, v\jiich are determined by the appro- 
priate nuclear charge and current densities. These form 
factors are in fact the Fourier-Bessel transforms into 
momentum-transfer space of the nuclear densities. Stated 
differently, the charge and current densities are the 
Fourier-Bessel transforms of the form factors into the 
co-ordinate space. In Section 11,3 we discuss the outline 
of the procedure involved in the reconstruction of the 
densities (from the available experimental data) in terms 
of their Fourier-Bessel parameterization. We next discuss 
in Section 11,4 and 11,5 the use of projected HFB intrinsic 
states for the calculation of the electroexcitation form 
factors. Although the projected HFB method has been 
extensively used in the past for the calculation of sub- 
shell occupation numbers, level energies^ static quadrupole 
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moments and reduced E2 transition probabilities, its use 
vis-^-vis the calculation of the electroexcitation form 
factors as well as transition densities has not been 
reported in the literature. 

11,2 Electroexcitation Form Factors for Nuclear Levels 

We shall be working in the perturbation theoretical 

1 2 

framework of scattering * , The electron- nucleus inter- 
action may be considered as an exchange of virtual photons 
and can be represented by the Feynman graphs given in 
Figure 11,1, Here Figure II, 1(a) shows one-photon exchange, 
where photon is ccupled to the electron by the coupling 

constant e and to the nucleus by the constant Ze, So the 

2 

expansion parameter of the series is Ze . This graph 

gives the first Born amplitude contribution. Graph of 
Figure 11,1 (b) gj-ves the second-order Born amplitude where 
two photons have been exchanged and the contribution is of 
the order of (Za) , The perturbation series is a power 
series in powers of (Za) for the amplitude. In general 
the cross section involves square of the amplitude. The 
lowest term in cross section is proportional to the square 
of graph of Figure 11,1 (a) and is known as the "fiirst Born 
approximation" of the cross section [having order (Za)^], 
Next higiest term contributing to cross-section is of 
order (Za) coming from the cross product of Figure II, 1(a) 
and Figure II, 1(b) which is called the second Born 




approximation. Square of Figure II, 1(b) and the other 
graphs contribute to higher power of ( Za) , In all these 
cases initial and final particles are considered free. 

In the 2p-1f shell nuclei, vvith Za 0,1 4-0, 21, the contri- 
bution of graphs other than the one-photon exchange ones 
are expected to be small. In what follov-/s, we discuss 
nuclear electroexcitation in the one-photon exchange 
appro ximation. 

Consider an electron with four-momentum 
incident on a nucleus of mass M, After exchange of virtual 
photon with four-momentum A = (q, iE) the electron gets 
scattered carrying four-momentum 

electrons vdth the energy of the order of a few hundred.. 
MeV, we can take and E 2 ^ k 2 . The different kine— 

matical relations satisfied by these quantities are 


I 

ro 

H 

q = (k^ + k 2 - 2 k;|k 2 Cos«)^^^ ; 

(II. 1) 


k^+M = k2 + E’ 

(11^2) 

where 

E’ = ( q^ + 

(II.3) 


Here M’ is the mass of the excited nucleus. 

Let the quantities (A^, i^g) characterize the electro- 
magnetic field of the passing electron. The quantities 
summarizing the electromagnetic properties of the nucleus 
are; the nuclear charge density eP(r ) normalized to 
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the nuclear current density, e3^(r) and the nuclear .. .. 

w 

magnetization density ep.^( r) ^ The interaction Hamiltonian 
which completely describes the electroexcitation process 
is given by 


Hj = He + Hj. . (11,6) 

where the Coulomb part is 

^c * , (11,6) 

and the transverse part is given by 

Hf = -e/[ J^Cr) ^g(r) + (I^C r) . 5x ^^(r) ]d®r, 

(11.7) 

The transverse interaction plays role' mainly for 
large angle scatterings. From now onwards we shall be 
concentrating here on the Coulomb part only. 

The Coulomb potential of electron is given by 

«e(r) = /[ Pg(r’)/lr - r’ l]d^r» , (II .8) 

So we have 

= e/[PNC^) PgCrM/ir- rM] d\ dV .(11*9) 

mmi -a ^ 

We write R » t/ir - r* | in terms of its Fourier transform; 

* 2ir^/ exp(iq,R) [ exp( -iq,!l’ )/R» ] d%*d^q. (II.IO) 
After integration w,r,t, R’ , we get 
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|r- rM~^ =(271:^)"'^/ exp[iq,(r- r’ ) ] q'’^d^q , ( 11,11 ) 

Now we make use of iHayleigh expansion, i.e^ multipole 
expansion of a plane wave 

exp(iq,r) = ^ i\2L+l) jL( Pj^( Cos 0^ * 

( 11 . 12 ) 

Using 

®12^ “ M ^ ^Lm(^)» (11.13) 

v\/e obtain 

exp(iq.r) = 471 

Now Coulomb interaction can be written as 

He = ^ ^( cj) » (11,15) 

where the nuclear Coulomb matrix element is given by 

-''f*|.j(r) (11.16) 

and X( q) is the Fourier transform of electron charge density, 
X(q) = /pg(r') exp(-i q,r» )d^r* , (11,17) 

Here we shall be making use of the plane wave solution , 
of Dirac equation for initial and final states of the electron, 
or in other words we shall be working in the framework of 
the plane wave Born approximation. 
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Initial and final states of the electron are 

— u( ) exp i(kj^,r - <i)t), = u(k 2 )exp i( k 2 *r— cot) 

(11.18) 

and the expression for transition charge density is 

Then expression (11,17) reduces to 
X( q) = ~e / u( k^ ) exp i(~k 2 .r’ +(i)t)exp i( k^ ,r-a)t) 

-* -* 3 

x exp (~i q, r' ) d r’ 

= ~e ( 2tx) ^ 6 ( k^ - k 2 - q) ( u( k 2 )Y ^ u( k^ ) ) . ( II ,20) 
Substituting this in (11,15) we have 

M 

The differential scattering cross section is obtained 
2 

by multiplying by 2n: tines the phase space factors. 

The phase space factor is given by 

(2tc)''^ d^k2 6(k^ + M - k2 + E') , (11.22) 

where normalizing volume (Q) has been taken to be equal to 

•* 2 

1 , Uberall has pointed out that the phase space factors 
are not affected by dynamac (nuclear) recoil effects involving 
free motion of nuclear centre of mass,. Making use of 
invariant normalization (m /E) for spinors u( k „) and dividing 

“ 1 f ^ 
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by the incident flux k^/E^, the cross section in the 
laboratory frame is given by 


dcr= (27iE^/k^) 6 ( k^+M-k 2 ,-E’ )( m^/E^ E^) IHJ ^ d^k 2 ( 2 tt)'^. 


(11.23) 


If we consider the scattering of the outgoing electron 


in the solid angle dQ, then 


/6( k^+M-k^-E’Ck^)) = k^ dQ H + dE’( k^dk2l"'’ 


^2 dQ[1 + (2k^/M)Sin^(©/2)]. 


(11.24) 


Since dE‘(k 2 )/dk 2 = ( dEVdq)( dq/dk 2 ) 

= (q/E*) (dq/dk 2 ) , 


(IIi25) 


where we have used the equations (11,3) and (11,1) and the 
approximation E’ ^ M, k^ = k^ . 

2 

The actual cross section involves a factor IHI 

c 

and hence the expression Iu(k 2 )YQ u(k^)l^'. For unpclarized 
beam this has to be summed over final states and averaged 
over initial states. We need to evaluate 


^ sf.Sg . 

/* 

IVlaking use of the following relation'' 


( 11^26) 


S "^6 »S2)%(k,S^) = [(K + , (11.27) 

, S2 



and the trace identities for even number of matrices 


•••• •^2n^ — ^1*^2 ^^^3*** 2n^ "" ^ 1 *^3 ^ 2^4 *'* * ^ 

+ (11.28) 


tr Y° Y ° = 4 , 


the expression (11,26) simplifies to 

Tr[(Y°()^2 

= ( E^.E^ + .k^ + ’^e^/®e • ^ 


Next we make an extension of the formalism v/ith a 
view to incorporate the possibility of nuclear electro- 
excitation, ’^e replace the nuclear matrix element 
by its operator definition 

“ i^/P°P(r) jL(qr) Yj^j^r) d^r (II, 30a) 

“ I ®k \M^^k^* (II, 30b) 

where P°^( r) » e^6( r - Tj^) , 

Let us denote initial and final states of the nucleus 
by fJ^M-^>and 1J^M^> respectively. Then the reduced matrix 
elements of fj^j^ are 

<JfM^lfLM = [Jfl'"’’ <Ji LMlJ^Mp< J^rif°PlIJ^ >'^ 


(11.31) 



(11.32) 


where we have made use of ?/igner--Eckart theorem. Here 
[a]-'' = (2a + . 


To obtain the differential cross Section (11,23), we 
collect the terms involved in it and substitute from (11,21), 
(11,24), (11,29) and (11,30), Starting with ground state • 
»J.M. > , the operators and H occurring in |H 1 can 
take it to final states iJ^M^> and IJ|rM|r> respectively. 

Then the differential cross section comes in the form of 
"density matrix" with indices (J^i^) and (J|M|,): 





= (4ii)^ (2e\2/N [ JfKJ^l) 

S ( J.M, ,LM !J.M.)( J.M, ,L 'M' U-^.) 
LM,L«M' ^ ^ 

X , (IIi33) 


where 

Cl = vye) y*,„,(9) \^t:) 4^ . (11.34) 

V(,(e) = + h.kj (11.35) 

and 

fL^(q) = <JfflfL(q)'*Ji> . (11.36) 


Finally, the differential cross section for electron 
scattering v/ith simultaneous nuclear transition ■* 
is obtained from the density matrix by putting J|. = 
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M|j » and summing over final projection and averaging 
over initial ones Carrying out these operations, we 

get 


Ji •* Jj 

/ dd ^ ^ . 

'^Coulomb 



%Mf M. 

1 f ’ f 


(11.37) 


This can be simplified by using orthogonality of 
Clebsch-Gordan coefficients 

2 (Jj Mj., LMIJ.M.)(JjM 3^,L»M« IJ Mf) = [J.]^ [L] 

\l' °mm« 

and 

f. \m^ * [ L] / 4it. 

M 

We finally get the expression 

= (27^6^2/^“^) V<,{e)IF(q)|2 . , 

(11.38) 

where the recoil factor is added for the sake of completeness 
and is given by 

^recoil ""1^^ +(2k^/M) Sin^ , (II,39a) 

and the form factor is given by 

IF(q) (4-Jf) ^ = (4ic/Z^)(2J^+l)“'‘ f, KJf Ji> 4 



(II. 39b) 



Z-O 


Further, in order to compare the form factors 
calculated by the PWBA with the experiments, the experi- 
mental data is plotted at the effective momentum transfer 
Qeff instead of the kinematic q, 

%ff = <? [1 

where is the radius of the equivalent uniform charge 
distribution of the ground state and is tho energy of 
the incident electrons. 

Sometimes the cross sections are given as a function of 
scattering angle. In these cases the momentum transfer is 
calculated using the relation 

q = 2(E^E2)^/^ Sin(0/2) . (II .41) 

For the evaluation of the single particle matrix 

elements we have employed tho oscillator vrave functions with 

the length parameter given by b = 1,01 A , The contre- 

2 

of-mass correction has been taken into account by multiplying 
2 2 4 

a factor exp( b q /4A) , The correction due to the finite 
size of proton has been incorporated by multiplying with 
the proton form factor f ( q) ; 

ir 

fp(q) = 1/(1 + q^ap/ 12 )^ , (11.42) 

where a = 0^84 fm; 



24 


II. 3 Nuclear Transition Charge Density 

The nuclear charge density operator P(r) appearing 
in equation (11,30a) can be expanded in multipoles; 

The charge density can be characterized in teims of 
its reduced matrix elements between the initial and final 
nuclear states of spins and (refs, 5-7) 

< I P °P( J.M^ > = 1 1 P °P( r) H J . > 

^ • 

(11.44) 

The nuclear electroexcitation form factor is related to 
the transition charge density by the relation 

FL(q) (Ji-Jf) =[(2J^+l)/(2Jj^+l)] / <Jf flP^P{|J.> 

o L 

X jL( qr) r^ dr , 

(11.45) 

The form factors are thus the Fourier-Bessel transfoicns 
into momentum-transfer space of the nuclear transition charge 
densities. In other words, the densities are the Fourier- 
Bessel transfoims of the form factors into coordinate space. 

In view of the relation ( 11,44} it is clear that the 
nuclear transition charge densities and the form factors are 
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equivalent quantities. However, the form factor acquires 
significance — - and can be considered more basic — in the 
case of the light nuclei (v/ith small Z) because the relia- 
bility of the PWBA in these nuclei permits one to relate 
directly to the experimentally measured cross sections 
(see Eq, 11,41), On the other hand, for nuclei with large 
Z, distorted wave Born approximation (DWBA) calculations 
become necessary since the Coulomb distortion of the electron 
wave function is expected to be significant. In the frame- 
work of the DVifBA. calculations, the experimentally measured 
cross section can be related to the nuclear transition charge 
densities in a form that cannot be presented through the 
form factors given by Eq, (11,44), The simple relationship 
between the cross section and the charge densities (via the 
form factor) is no longer present. In such cases the 
transition charge densities play a more basic role since they 
summarize the nuclear structure information independent of 
the precision of the PWBA, 

Fourier-Bessel Parametrization of Transition Charge Densities 

The basic assumption underlying the Fourier-Bessel 
7 

parametrization is that the transition charge density is 
restricted to a region r < and may be expanded into a 
Fourier-Bessel series with the spherical Bessel function jj^ 
of order L ; 
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00 

for r < 

0 

for r > R 


(11.46) 


Here the q ’s are given by the condition 





In the PWBA framework, the coefficients Ajj^ are 
easily expressed in terms of the form factors F^( qp^) 
[i.e, measured at the momentum transfer q^ ; 


A 


P’ 


2( 2Jj^+T ) 1 

t'2Jf+U Ro 


Fl(c^) .(11.47) 


Equation (11,47) has been derived by using the expre- 
ssion (11,45) in conjunction with the following relation: 

^c 

/ tj'i “ ■# tji+1 • 

( 11 ^ 48 ) 

In view of the non-availability of the form factors 

for momentumr- transfers q>q ,the Eq, (11,47) permits the 

max 

determination of only a certain number of the coefficients 
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+ 4 * 

Considering the 0 2 transition and taking 

= 1 0 fm, it is seen that the usual data (covering the 
range q = 0 - 2,5 fm ) permits the determination of seven 
coefficients, A (p,= 1,2,., .,7), corresponding to the zeros 
of occurring at qR^ = 5.76, 9,10, 12.4, 15.5, 18_i7, 

21,85 and 25,0. In the case of the O”*” -* 4’*’ transition, an 
identical range of the available experimental data facili- 
tates the evaluation of six coefficients corresponding to 
the zeros of 34 ( 9 ^^^^) at QR^^ = 8.18, 11,7, 15.0, 18.3, 21,5 
and 24,7, 

For the unknown large-q behaviour of the form factor 
for q > (<« 3 fm ) one usually makes the following choices 

dictated by general physical arguments; 

(a) The relation 

FL(q)<c q-"^ f p( q) (11,49) 

is considered as the upper limit for the large-q form factors. 
Here f p( Q) is the form factor of the proton ( eq, 11,42) , 

The constant c is adjusted by matching the observed form 
factor at the second or the third maximum, 

( b) The upper limit 

FL(q) < c fp(q) 

is also sometimes employed in the Fourler-Bessel analysis; 



Outline of the procedure employed for extracting the 
transition charge 'densities from tHe measurea form 
factors 

The procedure involves the use of, the PWBA as well 
as the DWBA methods in the following manner; 

( i) One first chooses a cut-off radius * 10 fm. This 

value is expected to be reasonable for the 2p-1f shell 

nuclei. The zeroth-order values of about 15-20 coefficients 

are calculated using phenomenological form factors (such 

8 

as the Tassie form factors ), 

( ii) One next determines the first 5-6 coefficients 
corresponding to those lying inside the experimentally 
observed momentuia-transf er range by invoking the KJfBA 
framework, 

( iii) The next step involves the use of_the transition 
charge densities obtained in the step (ii) to transform 
the experimental cross sections to "experimental” form 
factors F- • One employs the relation 

Li p 0Xp 

^L.exp “ '^L.Calc. ^ 

( 11 . 50 ) 

where ( is obtained via the DWBA framework using 

the Fourier-Bessel expansion of P ( r) obtained in the step 

(ii), Eq, (11,50) is easily derived by considering that 

the value of { is simply the square of 

and by assuming that the fractional difference between the 

experimental and the calculated PWBA cross sections ' 
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is the same as the fractional difference between the 

experimental and the calculated DVffiA cross sections, 

• 

( iv) The "experimental" form factors are next employed 
for carrying out an optimization of the complete set 
of expansion coefficients involved in the step (i). 

This is again performed in the PWBA framework. One also 
considers here the variation of the cut-off radius R 

c 

as well as the effects due to different choices for the 
large-q behaviour of the form factor. The optimized 
set of parameters is employed in repeating the procedure 
from step ( ii) onwards. 

The influences of the cut-off radius and of the 
large-q behaviour has been examined by Neuhaasen by 

considering, as an illustrative example, the case of the 

+ 64 

2^ state of Zn, 

11,4 The Pro jected Hartree-Fock— Boqoliubov Method 

In this section we discuss the calculation of the 
nuclear matrix element appearing in eq, (II, 39b) in terms 
of the states IJ^> and |J£> obtained via an explicit 
projection from the HFB intrinsic states. In this connection 
Vife first discuss the basic elements of the HFB theory, 

II,4,1 The HFB Method 

The many-body Hamiltonian of the nucleus can be 


written as 
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H = S <o<UIa >a;a^. 1 <ap,v^ I T 6 >a;a; a^a^ (11.51) 

with summation indices running over the orbits constituting 

the valence space. Here ^ and a^ are particle creation 

and annihilation operators, respectively. The single 

particle energies, <alela>, include the contribution due 

40 

to the spherical single-particle field of the Ca core. 
Further, <a§ 1 lY6>is the antisymmetrized two-body matrix 
element. In the HFB procedure, the first step is to go to 
the quasiparticle bases by applying unitary transformation 
on the particle bases and finally assuming that the inter- 

Q 

action between quasiparticles is relatively weak , 

We can write the Hamiltonian in terms of quasi- 
particles as 

H = + V + Vlnt> 

where is the energy of the quasiparticle vacuum, H^p 
describes the elementary quasiparticle excitations and 
Hqp^int “the (relatively) vfeak interaction between quasi- 
particles , 

The Hamiltonian (11,52) has many symmetries associated 
with it and correspondingly we have conserved quantum 
numbers. But the part ( E^ + H^p) of Hamiltonian may not, 
in general, preserve all the symmetries of H, In HFB one 
imposes the constraints by introducing Lagrange multipliers 



in such a way that atleast the obseorvables corresponding 
to the violated symmetries have the desired expectation 
values. More specifically, in the HFB framework one 
includes pairing correlations by introducing number non- 
conserving wave functions. The Hamiltonian H is replaced 
by 

h' = H , CII;53) 

where the Lagrange multipliers and are chosen in such 
a way that the number operators 

i* i* 

=s 2 a a and N = S a .. a 
^ a ocTc ait V a av a v 

have the expectation values 

< $0 I ^7: * ^ 0 ^ ^ and ^ I Nj, > = A-Z , 

(11.54) 

The general Bogoliubov transformation employed in 
the HFB method is 

^ • ( 11 - 55 ) 

So each quasiparticle is a linear combination of all 
particle creation and annhilation operators. U and V are 
NxN complex matrices for N single particle basis states. Then 
2 Nx 2 N linear transfer nation 



q 

q 


r 


U V 


V 


U 


\ 

f 4 . ^ 


t 


a 


a 

J 



(11.56) 


is required to be unitary. The unitarity of the transfoi>- 
mation matrix 


i.e, MM^ = 


and 




< 

u 

N 

V • 



M 

S3 



$ 

(11.57) 



V 

u . 



u\\ 

I, 

implies the relations 


t 

f 


f 

* 

( 11458) 

uu + 

w 

33 

u u 

+ V V = I 




+ 

~ * 

(ii;59) 

uv + 

vu 

33 

uv 

+ V U = 0 


where , U and U stand for the adjoint, transpose and 
complex conjugate of the matrix U, These relations can 
also be satisfied by requiring that the particles and quasi- 
particles obey the feimion anticommutation relations 

C4 » t ^ 0 

(11.60) 






(11.61) 


Since M*"^ *= , we can invert (11,55) to express 

particle operators in teims of the quasiparticle operators 


a‘^=2i(u„ qt+^4 qo) 

a p ^ pa p pa ^ p ^ 


(11.62) 
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The quasiparticle vacuum is defined through the 
relation 

q = 0 for all a (11,63) 

cc ° 

which has a solution 

]§> = N 2q }0> 

^ a cc 

where N is a normalization constant and ! 0 > is the particle 
vacuum. 


Derivation of the HFB Equations 

In the derivation of HFB equations we shall be making 


use of the expectation values of the operators a^ and 

a c 
a 

by 


a^a^. We define density matrix P and the pairing tensor t 


'’cp = < *0 • I *0 > 


(11.64) 

(11.65) 


Using equations (11,61) and (11,62) we can write 


‘’ap = 




i.e, P 


S S V V* 6^. = S V.„V* 
y 5 yp 6a 6 6p 6a 

t 

V V , 


(v'^V) 


ap 


( 11 . 66 ) 


and 





< *o'[I [ S ( 'V"+ 

= ? f “ep'^ya ®«Y = i ''6ccUsp= (v"u)„p 

or t = V^^U . (11.67) 

Using equation (11,59) we can write 

t = (u'’'v)'^= -(V U*)’*' = ~U V* = -t . (11.68) 

It follows from equations (11,66) and (11,68) that P 
is hermitian and t is antisymmetric; 

As a next step, we are interested in the expansion 
of operators appearing in H’ in terms of normal ordered 
products. Normal product of an operator (made up of 
products of annihilation and creation operators) is obtained 
by expressing the particle operators in terms of the quasi- 
particle operators and then rearranging them in such a way 
that creation operators are to the left of the annihilation 
operators, A minus sign is included for an odd number of 
operator permutations. The contraction <0> of an operator 
is defined as the expectation value of 0 with respect to 
quasipartacle vacuum. 

The expansion of operators in terms of normal products 
is given by Wick’s theorem which states that any product 
of operators 0 is equivalent to the sum of normal products: 0: 
which contains all possible contractions. The Hamiltonian H’ 
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involves operator products of the type p and 
From Wick’s theorem v;e can expand these operators as 


*®a ^p* ( Ii*69) 

and 



<^a 

ay > 

<a'^p a5>+<a^ 

4 > 




< a 

Y 



+ 


• ^ a ® 6" ” 

^ > 




• a"^ 

. ap 



4- 

<^a 

ap^> : 

®Y ^6° + < Sy ag > - 




5 

B- 


+ :a + a+aYag; {[1.70) 


Substituting (11,69) and (11,70) into the equation 
(11,53) and using definitions (11,64) and (11,65) we can write 
H' as 


H’ = ^ ^2 ^4 (11*71) 

I 

v\here H has n uncontracted operators: 
n 

= tr [ (e- (l/2)r)P+ (i/2)At'^], (11,72) 

^2 ^ ap “\-N^)ap ^p* 

(11.73) 
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( lli74) 

Here the HF Hamiltonian h, the HF potential r and 
the pair potential A are given by 

h = e+ r , cay I ip6> (11.75) 

and A^p = > ■''YS (11.76) 


By invoking the properties of < f V^! >: 

< a p ! Iy 6 > = - < pa I I Y 5 > j < ay I I p 6 > = <p 6 lv^[aY>, 

and using the hermiticity of P and antisymmetry of t, it 
follows that h and r are hermitian and A is antisymmetric. 

Since the expectation values of all normal order 
products involved in H’ vanish, we have 


< §0 IH J§o> = tr [( -Xy + ir) P + i A t^ ] 

( 11,77) 

I 

The H2 part of the Hamiltonian describes quasiparticle exci- 

j 

tations and the remaining H^ contains quesiparticle inter- 
actions . 

I 

In HFB approximation, H should acquire the form 


H 


,H^+ 2 Eq^q + H .. 

o „ a a a qp-int 


(11.78) 


This can be achieved if we assume that the ’’independent 


quasiparticle” Hamiltonian has the form 
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= S E q’’' q 
^ a a a a 


(IIi79) 


I + 

The commutator [ H„, q ] is given by 

^ tx 


[ ^2’ ^a~ 


®p^ • (II.SO) 


This commutator evaluated by using equations (11,73) 


yields 


[ q^]= 2[h' U + A„ V lat 

■- 2' ^a-* pY PY “Y PY “Y ^ P 


+ 2 r 
p Y 


%■ 


^pY ^aY^^P * ( ) 


where h = e ~ X - X 

% V 


On comparison of the equations (11,80) and (11,81) we get 
the general HFB equations: 


= Ei 


(11.82) 


where = (U^^» U^2» •••••» ’^iN^ • 


. ime-rev i 


timetry 


We have assumed time-reversal symmetry in our 
calculations in the 2p-1f shell nuclei. All the basis 
states in the configuration space spanned by the orbits 
1 ^ 7 / 2 * ^^3/2* ^^1/2 ^%/2 divided into two sets. 

The first set containing states I k > is characterized by an 



even value for (m + 1/2), The set of time-reversed states 
I k > = T Ik >has (m + 1/2) = odd integer. The phase conven- 
tion is; T I nljm > = (~1 )^’''^“"^In Ij-m >; 

The states lk> are taken to be: 


(1f7/2 7/2)» ^‘•^7/2 3/2» ^^5/2 3/2»^^3/2 3/2^ 


^^^7/2-1/2» ^%/2-1/2* ^P3/2-1/2» ^Pi/2-1/2^ 
and (1^'7/2-5/2’ ''^5/2-^5/2^ • 

For qua si parti ole transformations 



^ (U „ 
p ap 


(U. 







) 

) 


(11.83) 

(11.84) 


we immediately notice that P cannot connect 1 k>and I k > 
states and t has non-zero contribution only for elements 
connecting 1 k> and lk> so we can partition P and t as 



(11^85) 


Since the interaction V conserves the magnetic 
projections the potentials (r,A) can also be partitioned 
as ^ 



f 






0 


0 

h = 


> 

II 


0 

h^ 


A Q , 

2 ^ 



«« 




h 


( 11 . 86 ) 



Time-reversal symmetry is imposed by requiring 

= T qt , (11.87) 

a ^ 

so that , v^p = - V*p (11.88) 

The quasiparticle vacuum is now invariant under time- reversal* 


Canonical Representation for the Time-reversally. Symmetric 

^ 


Next we are interested in the simultaneous diagonali— 
1 0 

zation of P and t. Creation and annihilation operators 
can be treated simultaneously by defining 



(11.89) 


Here A"^ contains quasipapticle operators and has only 
particle operators. Then the general quasiparticle trans- 
formation is given by 


a'*' = M . ( 11.90) 

We can now define the generalized quasiparticle density 
matrix Q by 


Fran the commutation relations of quasiparticles (11,61) 

we have 




( 11 . 92 ) 
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Substitution of (11.89) and (11,92) into (11.91) yields 


Q 



(11.93) 


where ’0’ here is NxN null matrix and '1' stands for NxN 
unit matrix. 

Similarly we define the generalized particle density 
matrix R by 



< *0 ' 4 Ba 



C 11.94) 


Substitution of (11.89), (11.64) and (11,65) into (11,94) 


leads to 


R 



t 


1-pJ 


(11.95) 


By looking at the inverted relation corresponding to 
(11,90) and comparing ( 11,91 ) and (11,94) we find that the 
generalized particle and quasiparticle matrices are related 
by a unitary transformation 


R = M'^QM , (11,96) 

Since m'^M = 1 and = Q, it follows that 

R^ = (m'*’ QA)^ = m’^Q(M M‘'')Qf4 = QM . (11,97) 

2 

The relation R = R yields from (11,95) 


+ t t^ 


P 


(11.98) 
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and Pt + t- tP = t orPt = tP 


(11.99) 


Thus ^ -j 't-j = This ensures us that 9^ and t^ 

can be diagonalized si:niultaneouslyi. 

We can now have a basis spanned by[Ik^> > Ik 2 >, .... 

I k^ > , f lc 2 > ] such that P and t are in the 

canonical forms 



where the basis states are taken to be j k >= bj^ T 0 > 


Let and t^ have the eigenvalues Py. and t 

KK 

respectively then from (11,98) v;e obtain 

It I (11.102) 

kl ^ ^ 

putting Py = V^, ^k ” 


It I*= . (11.103) 

kk ^ 


To satisfy the equations (11,102) and (11.103), the 
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quasiparticle vacuum can be expressed as 


I > = 


I (Uk + Vk b+ b2) ! 0 > (11.104) 


Approximations employed in the present set of calculations 
We recall the special quasiparticle transformations 




q; = U* b t _ v; b^ 


and 


<’2 = + 


(11.105) 


(II;106) 


After inversion of this transformation and substitution 
into two-operator part of H’ we can expand as 

4 " ^11 4c (11.107) 


with - ^kt \ 


^1*1 


H 


20 


kki ^^bO^kk’ ^k» ^k ^ (11.108) 


'20 ^kk’ ^k> 


20" kk’ 


In order to realize the free quasiparticle approximation, 
as we have seen earlier, the Hamiltonian (11,107) must satisfy, 

I 1 

the condition H 2 Q = 0, and H-i-i must be diagonal. The condition 
(H«n) = 0 leads to the BCS equations and the requirements 

20" k"^ 

(^ 20 ) - = Of ® ^ ^ ^ gives HF - like equations 

klc 


t 


provided h,. , = 0 and A _, = 0 (with k k' ) . We have 

kk 

considered here the latter approximation, viz. A ,= 0, The 

klc 



coefficients ^ appearing in eq. (11,101) result frcan a 

diagonalization of the HF - like potential h’ (v;hich includes 

2 

the relevant density Vj^) ; 


(IIi109) 


The occupation probabilitiies are obtained from BCS 
equations; 

_ 1 S 


..1 


i k' ^kkfV^l k» k > Uj^, V^, . (IlillO) 


Let be the eigenvalues of (11^109). Then the 

f 

condition (H 2 q) „ = 0 leads to 


kk 


2«k - A _ (U? - V?) = 


kk 


k 


( 11 . 111 ) 


where 


2 . 2 

= 1 


( 11 . 112 ) 


The set of equations (11,111) and (IIi112) can be solved 
easily by making substitution 


V 


= Cos 9 and Uj^. = 55 n 9 


(11.113) 


From (11,111) we have 


Sin2 9 = - A Cos 2 9 

or tan 2 9 = “ ^/®k ’ 

or Cos 29 = 


(11.114) 

(11.115) 



This equation together vdth (11,112) yields 


= 5[1+ V(®k (11.116) 

and 

(11.117) 

Ick 

Substitution of equations (11,116) and (11,117) into (11,120) 
gives 

’ 2 k‘ < ^ 

(II;118) 

The condition 2 ^ ^k ^(v) ” ^ 71 ( 1 ^) 

I [1 - ei/(e2+Ag^)i/2] = (11.119) 

%{V) 

To discuss the iterative procedure used in the 
calculations vve define the quantities 

1 0 

PEP-t ’ ^22 ’•***^0T0^ ^ 

( 11 . 120 ) 

Q[ii»(A,Y » ^22 » ••*•^0 To^^ = (i/[(®k 

( 11 , 121 ) 

Here ^ and ©j^ = ^y. ~ P , (11,122) 



Summing both sides of eq, (11.117) over k, v;e have 

^^■n/v = 10.0- P + Q (II;123) 

We thus obtain a convenient expression for p. t 

p = (P - 10.0 + \/y)/^ . (IIi124) 

We have follo'.ved the following iteration procedure 
in our works 

a. We first carry out the usual HF calculation. This 

provides an "educated” initial guess for the set of 
values of (k = 1,10), p as well as the expansion 

coefficients defining the orbits I k> (k = 1,10), 

b. We next assign a constant non-zero value (say 0,5).. 

to all The equations (11,120), (11,121) and 

(11,124) are iterated until one obtains a "self- 
consistent" value for p , 

c. This value of p is next employed to obtain a set of 

self-consistent values of A ’ s by using repeatedly 

klK 

the equation (11,118), 

d. The A values resulting from the preceding step 

kk 

are next employed to obtain an improved self-consistent 
value of p by again using the equations (11,120) and 
(11.121) and (11,124). 



e. One obtains finally a mutually consistent set of 

values of {J. as well as A 's by repeating the 

kk 

preceding steps (c) and ( d) , 

f. We next solve the equation (11,109) to obtain a 

new set of and iV, The procedure just discu- 

ssed is then repeated from the step c onwards until 
successive iteration cycles produce nearly the s_ame, 
set of numbers for various parameters. 

In the HF-BCS procedure one simply terminates the 
calculations at the step e. 


1^*4, 2 The Technique of Angular Momentum Projection 

We discuss in this section the general technique^^ 
of projecting the intrinsic deformed wave function (of 
the HF or the HFB type) on the subsoace of angular momentum 
J, The intrinsic state involving axially symmetric deformed 
single-particle orbitals is characterized by the eigenvalue 
of Jj, operator. We denote the intrinsic state by * ^ 

(K = ^ expand it in terms of states of good 

angular momentum; 


i > = 

K 


S 

J 




> 


( 11,125) 


Applying the rotation operator 


R( Q ) = e 


-iaj 


-iGJ, 


-iYJ. 


(11,126) 
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we obtain 


I 1. I > 


JM J n-MC 


(11.127) 


where are the matrix elements of the rotation operator 
in the basis IlPv> • Q stands for Eulerian angles a, Q 
and Y • 

.J 


Multiplying (11,128) by and integrating over Q , 


we have 


^i> = 


2 J + 1 


/dQD^j^(Q) R(fi) > ( 11 , 128 ) 


871 a 


\Miere we have used 


2 TE tt 27E 

/ dQ = / da/ dG SinG / d Y 
O 0 o 


( 11.129) 


and the orthogonality of the D-ma trices: 


/ d Q ( Q ) Q ^ " 2J+1 ^JJ ’ ^MM' (11.130) 


8 




In view of the equation (11,128), we can say that 
the state of angular momentum J v;ith projection M in the 
laboratory frame is projected out from the deformed intrinsic 
HF/HFB state 1 ^ possessing projection K of J on the 
Z-axis of the intrinsic frame. 

Equation (11,128) can be written as 

I 5 ^K> (11.131) 

wh ere 
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” 8?^ •'■dQD^(Q) R(Q) 


( 11 . 132 ) 


.J 


The projection operator satisfies the following 
properti es; 


" ^I<K 

(ii) (PkK^ ~ ^I<K 

(iii) [H,P^]= 0 


( 11 . 133 ) 

( 11 . 134 ) 

( 11 . 135 ) 


Making use of these properties of the projection 
operator along with (11.131), the energy of the projected 


state 


Ej = 


< IHI 

T“ — 

^K> 


( 11 . 136 ) 


can be expressed as 


% ~ 


<5 ! K P'^ F i > 

K ' ^KK ' K 


Using the expression for P^ and R( Q) given by 
equations (11,132) and (11,126), respectively, v./e have 


( 11 . 137 ) 


n j 


/ d© Sin© djQX©)< IH e 


I 


7 d© "sin© "d^Tek ^ ! e''^^y i > 


( 11 . 138 ) 


K 


where we have used the. property that exp(~iaj^) and 
exp(-*iYJ 2 ) operating on <®j^I and I states produce 
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exp( -iocK) and (~iYK), respectively, 0) are related 

to by 


D^XO) 


dJ^Xe) 


(11.139) 


11,5 Electroexcitation Form Factors and Transition Charge 
TSensiti^ iri of Frojected ¥avGfunc'€ions 

The axially synmetric intrinsic deformed HFB state 

1 2 

with K=0 can be written as 


I «o> = 


(u? + v? bt bt ) I 0 >, 
im ^ 1 1 im ! * 


( 11.140) 


where 




2 

j 


,ni 


a^l and b'^_. = 
im 


2 f -1 'l 3”^ a 

3 ^ ji * 

(11.141) 


The index i distinguishes the states with same m 
value and j labels the spherical single particle orbitals 

2p^y2» ^^3/2' ^^7/2* 


We can rewrite 1 in the form 


J 5q> — ^ 10^ ( 11.142) 


with 




m 




= E C 


a 




oy, 

V 


V, 


6 

ma ra 


Ui 


a 


.-“s 


(11.143) 
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Here “ stands for quantum numbers (Jq., %) and N 
is a normaliaation constant. 

The state .vith angular momentum J is given by 

= 2J^/D^^(a) R(0)i *^>dQ 


KK K ~ QT^- KK 


K' 


(11.144) 


Making use of equation (11,127) for R(Q) , 


" N/ dQD^Q exp (S^ F^^(Q) a+ a^") I 0 > 


ap ap 


a -p' 


(11.145) 


where 


F Q = 2 

aP m 


Ja 


.J'rt 

(11.146) 


J J 

Expression for overlap given by (apart from 


normalizations N) 


< > 
o o 


/Sin© d© d'J (©) < Olexp(i S^f* a„a ) 

Q 00^ ' 2 ap ap p 

X exp('2 P’ ^ ® p I ) 1 0> ( H.147) 

exp(^ Tr [ln(l + M(9))] 

[det ( 1 + M(9))]^'^^ (11.148) 


v\^ere 


M(©) F(©) f+ 


(11.149) 
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The general expression for the Matrix elements of 


13 


an irreducible tensor is given by 




J L J’l 
0 0 0 


X i: 


J L J» 
0 V- \i- 


^ d2D^(Q)<5^[R( Q) I0^> 


.L 

p ‘i-o' 


( II;150) 


14 


Any matrix element can be derived from equation 
(11,148) by differentiating with respect to f* and F, 

The expression 


< 5 JR(Q) o'P Itp e^pC 2 y^ 6^6 


(11.151) 


where 




2^<pfT^ lq>a+ aq 


then leads to the following final expression for the matrix 
elements of •* 


J- J. J. J.-1/2it/2 

< 1? ^ f f I ^ = [n ^ n / 2 

o LM=0 o ^ o P 


Jf L J^i 

-\i li 0 


d^^(G)n(r 


( 11,152) 

where the density matrix is given by 

P(e) = IQGZ&i '. (11.153) 
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and the normalizationsare given by 


n = / n( 0) ( 0) Sin0 d0| n( 0) 

o 


= [det(l+ M(0))3''^^ 


(11.154) 


The matrix elements of the Bessel function in the 

-f 

oscillator basis are given by 

<n’l’l lnl> = exp( -y)( ( n'-l) l( n-1 ) 

X (r(n +1 + '5) ^ ( n+l+ 2^)"''^^ 


n’-l 

X 2 

m'=0 


(- 1 ) 


in +m 


^ * 1 .. - 

m=0 m imi( n -m’-l )i(n-m-1 ) » 

r( ^ 1 ' +l+ 2 m' -.- 2 m+L+ 3 ) ) 

X . 

r(m + 1' + ^)r(m+H- 

O 

X F(^L-l»-l-2m’-2m ?L + y) 


( 11.155) 

2 

where y = (bq/2) ' and ? is the confluent hypergeometric 
function. 

The expression for the reduced transition’ probability 
for electric 2 -pole transition, B(EL, 0 ** 2 ), and the 

static quadrupole moments for the state, Q , , are given 

L 


B(EL, 0-^- th = :j^l<S^n^llS°>|2 


and 
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Qj^H- =<^ L ^ Qq l^L ^ (11.157) 

wh ere 

' (11.158) 

The expression for the reduced matrix element 
appearing in equation (11,156) can be obtained by replacing 
jL( by in the equation (11,152). 

The steps involved in the calculations are as 
follows; First the wave functions ) are obtained by 

carrying out the HFB (or HF— BCS) calculations. The ampli- 
tudes, (U^, V^), and the expansion coefficients, 
resulting from self-consistent calculations are used for 
computing the f matrices defined by equation (11,143), One 
next evaluates the matrices F, (1+M) andM(1+M) ^ appearing 
in equations (11.146), (11,149) and (11.153), respectively, 
for twenty values of the Gaussian quadrature points in 
the range (0,n:/2). Finally one computes the matrix element 
(11,152) in terms of these quantities. 

The transition charge density ^^(i') is the reduced 
matrix element of 

Pj^(r) =<J^nP°Pnji> (11.159) 

Making use of the equations (11,44) and ( 11,45) and 
the projected HFB wave functions (11,145), we obtain the 
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following expression for the transition charge density; 


J. J* 


Jf) = (n ^ n ( 2J.+1 ^ 


Ji L 
-p p 0 


a® 




X [M( 0) (1 + M( 0) ) ”"’ ] Sin0 d0 


(11.160) 


where Rf^]_(r) is the radial part of the harmonic 
states I nl > , 


oscillator 
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CHAPTm III 



Introducti on 

Recent inelastic electron scattering experiments 

-f “‘r 

have provided valuable data on the 0 -* 2^ as well as 

0^ •* 4^ transitions in a number of doubly even 2p-1f 

shell nuclei , Wi-diin a measured momentum- transfer 

range upto 3 fm ^ the observed C2 form factors are 

characterized by two distinct maxima appearing at 

and 1 ,6 fmT^ , The maxima of the C4 form factors occur at 

-1 

q «« 1 ,2 and 2,2 fm , The relative magnitude of the fom 
factors at the two maxima is expected to provide a sensitive 
test of the nuclear microscopic models involved in the 
calculati ons , 

In this Chapter v/e demonstrate that the yrast wave- 

TC -}- + 

functions with J = 0,2 projected from, the HFB 

intrinsic states resulting from realistic effective inter- 
actions operating in the 2p-1f shell permit a fairly satis- 
factory interpretation of the available data on the 2^ 
form factors in some doubly even Ni and Zn isotopes (with 
58 < A < 68) in the second half of the 2p-1f shell. In 
this context we also discuss the microscopic prescription for 


calculating the transition charge densities in terms of 
the self-consistent projected HFB wavefunctions , 

Recently Mooy and Glaudemans^ have calculated the 
form factors as well as the transition charge densities for 
the 0 -* 2^ transitions in the nuclei In terms 

of shell model wave functions resulting from semi-empirical 
effective interactions operating in the restricted valence 
spaces involving the 

( 2p2/2*' ^^5/2* ^^1/2^^"*^^^ configurations. The calculations 
involved an independent variation of the effective charges 
for protons and neutrons so as to fit the B( E2; 0^ -* 2^) 
values as well as the ratio of the form factor at the first 
and the second maxine for each isotope. The procedure 
yielded significantly mass-dependent effective chargosj the 
values for the proton effective charge, e^, turned out to 
be (2.60, 3,30, 1.85) for the nuclei respectively. 

The available data on the reduced E2 transition proba- 
bilities in the Zn isotopes suggests enhanced rotational 

collectivity in these isotopes compared to that in the Ni 
8 9 

isotopes * , It is thus desirable to have an extended shell 

model description of these isotopes involving more than just 

one (lf7/2) hole. In view of the near-intractability of 

such a description - the dimens iomli ties of _ the matrices 

6 

involved in the earlier shell model calculations with 
smaller number of valence particles are already above 1000 — 
the available form factor data in doubly even Zn isotopes 


has not been analyzed in a microscopic framework thus far. 

In the present work we show that the use of J = 0^f2^ 
states projected from the HFB wave functions generated in 

2p-1f space leads to a fairly satisfactory unified 

interpretation of the available IFI^ ( 2|) in the nuclei 

58,60,62... , 64,66,68., . ^ 

and Zn in terms of reasonable value of 

nearly mass-independent effective charges. This set of 

effective charges is also shown to yield excellent quantitative 

agreement between the calculated and the observed values of 

the reduced E2 transition probabilities in these nuclei. 


In recent years considerable attention has been 
devoted to the coordinate-space reconstruction of the transi- 
tion charge densities associated with the electro excitation 

of nuclear levels in the framework of the Fourier-Bessel 

1 0 

method of Dreher et al . ; these analyses have been carried 

1 2 
out by Heisenberg for the Kfi isotopes and by Neuhausen for 

the Zn isotopes. The empirical transition densities acquire 

significance vis-a-vis the test of model wave functions since 

they contain structural infornation unaffected by the 

accuracy of the usual plane-vrave Bom approximation ( PWBA) 

for large momentum transfers. We have calculated here the 

transition charge densities for various Ni. and Zn isotopes 

in terms of the projected. HFB v^ave functions. By combining 

the model space transition densities resulting from the 

self-consistent wave functions with each of the two existing 

models for the core-polarization transition density we have 



attempted to examine whether the empirical values point 
towards a preferred choice. 

In Section III, 2 v;e present some details of the 

calculational framework. In Section III ,3 we discuss some | 

general qualitative aspects of the electroexcitation form 

factors associated with the 2^ levels in the 2p-1f shell 

nuclei,. In Section III, 4 we discuss the form factors, 

transition densities as well, as the E2 transition strengths 
+ + 

for the 0 -* 2| transition in some doubly even Ni and Zn 

isotopes. Finally Section III,5 contains some concludir^ 

remarks, i 


III, 2 Calculational Details 
111,2,1 Effective Interaction 

In the calculations presented here we have employed 
as an effective interaction for the ( 1 ^ 7 / 2 * ^P 3 / 2 '^Pl/ 2 *^ ^ 5 / 2 ^ 

space a slightly modified version of the Kuo-Brown (KB) 

11 12 
interaction suggested by fAcGrory, Wildenthal and Halbert , 

The KB interaction consists of the bare matrix elements 

derived from the phenomenological Hamada-Johnston potential 

and the renormalization due to the thre^particle one— hole 

(3p-1h) excitation in the Ca core. Me Grory, Wildenthal 

and Halbert^^ found it necessary to modify the T=1 part of 


the KB interaction for optimizing the agreement between 

42—50 

the shell model and the observed spectra for the O 


The modifications consisted of lowering the 


nuc 1 ei , 


matrix 


centre-of-gravlty of the <( Vl( 1 > J T=i 

elements by nearly 0,3 MeV, as well as raising the centre- 

of -gravity of the i )IVI(if^^2)( ^ 

matrix elanents by nearly the same amount. This modified 
version of the KB interaction is hereafter labelled as the 
MWH interaction. In Appendix A we have tabulated the two- 
body matrix elements of the MWH interaction. 

The single-particle energies characterizing the one- 
body part of the Hamiltonian employed in the present calcu- 
lation are (in MeV) ; ^( 1 ^ 7 / 2 ^ “ ^ ^ 

111,2,2 The Choice of Intrinsic States 

-}- * 4 " 

For calculating 0 2,j form factors in various 

2p-1 f shell nuclei, we have obtained the axially symmetric 
intrinsic states by following three different methods 
the HFB method for the nuclei ' Ni, the deformed HF-BCS 
method for the nuclei ^^i and ^^Zn and the usual deformed 
HF method for the nuclei The HFB and HF-BCS 

methods are discussed in Section 11,4, Pairing correlations 
between only the like particles are allowed in our HFB and 
HF-BCS calculations. The choice of the intrinsic state for 
a particular nucleus has been dictated by the following 
considerations. It is seen that the inclusion of pairing 
correlations (in the HFB framework) results in nearly spherical 
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intrinsic states in nuclei having either N=28 or N >34, 

This is related to the fact that the usual KFB method does 

not permit a treatment of the deformation and pairing 

degrees of freedom on the same footing because of the 

subshell closure at hl=28 and the approaching shell closure 

at rfc40? it tends to emphasize the latter degree of 

freedom. The HFB intrinsic states for the nuclei with 

N=28 or N > 34 therefore yield vanishingly snail ( < 10^) 

amplitudes for the yrast states with J > 2, These HFS 

intrinsic states are inappropriate vis-'h-vis the present 

calculation which reguires, as an essential prerequisite, 

a reasonably satisfactory simultaneous description of the 
TC -f- 4“ 4* 

J 0 , 2 and 4 states. The HF-.3CS intrinsic states 
for tho nuclei ^^Ni and ^"^Zn, and the HF intrinsic states 
for the nuclei contain = 2"^, 4"^ states with 

sizable ( > 10 %) amplitudes. 


ITI.3 Qualitative Features Associated with the Form Factors 
for the -* 2^ Transitions in'the'2p-lf Shell Nuclei 

In view of the relation (11,152) we expect the matrix 
elements < f [j 2 (qr)lf > and <flj 2 (qr)lp > to play important 
roles vis-li-vis the observed squared form factors. The gross 

qualitative features of the appearance^ of two maxima at 

—1 —I 

q 0,7 fm and q 1 ,7 fm can be understood in teiros'of 
the q-dependence of these matrix elements. We assume here 
that the effective charges do not have significant effect 
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on the q-dependence of these matrix elements. Figure III.1 
shows that the matrix element I<f (j 2 (qr) I f> I possesses 
maxima at q 0,8 and q 1 ,9 fm“^ , and the maxima 

of the matrix elements I <flj 2 (qr)?p> I occur at q ~ 0,7 fm-”* 
and q <'» 1 ,6 fm , The zeros of the matrix elements 

<f lj 2 ^ <fl j 2 (c[r)ip> are at q 1 ,6 fm"^ and 

-1 

q 1 ,2 fm , respectively. Since the expression for form 
factors involves a linear combination of the matrix elements 
<flj 2 (qr)lf> and <f}j 2 (qr)l p> , we expect the maxima 
of the form factor to occur in the ranges q = (0,7, 0,8) fm“^ 

— “I 

and q =s (1,6, 1,9) fm , Similarly, we also expect the form 
factors to display a zero in the range q = (1,2, 1,6) fm~^ , 
These qualitative features of the form factors are consistent 
with the characteristics of the observed form factors in 
various 2p-1f shell nuclei; 


III, 4 The o'*” ** 2^ Transition in the Nuclei 


Ni 

and V ' Zn 
111,4,1 Form Factors 

\'‘!e first discuss here the Ni isotopes. The squared 
form factors I F( q) I ( 0 2^ ) in these isotopes are calcu- 


J 


lated and compared with the experiments in Figure III ,2, 

The form factors have been computed vdth the effective 
charges e =1,8 for protons and e,, = 0,8 for neutrons. 

As discussed later, it is seen that this set of effective 

I I 

charges is quite consistent with the available B(E2j 0 *♦2^) 





FIG. II I. 2 EXPERIMENTAL AND CALCULATED SQUARED FORM FACTORS 
If]^ for THE o'*’ + 2 * TRANSITIONS IN THE NUCLEI 
58,60,62jj^^ THE SOLID CURVES SHOW THE RESULTS 
OBTAINED WITH THE SELF-CONSISTENT WAVE FUNCTIONS 
THE BROKEN CURVES DISPLAY THE RESULTS OBTAINED 
IN THE RECENT SHELL MODEL CALCULATION INVOLVING 
RESTRICTED COLTIGURATION MIXING. 
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values in these isotopes. 


Overall, it is seen that the agreement between the 
caxculated and the observed form factors is remarkably 
good throughout the momentum-transfer range q = 0,5-2, 8 fm"^ , 
The degree of quantitative agreement in the nuclei 
is quite satisfying, particularly in view of the fact that 
the effective charges were not fine-tuned for each isotope^ 

A significant discrepancy occurs in the case of the 
nucleus Ni where the present calculation’ underestimates 
the form factors for the range q = 0,8 - 1,2 fm“^ by about 
30 percent. This may be due to the non-inclusion of the 

1g9/2 orbit in the valence space. As pointed out by 

13 14 

Delphini and Glaudemans as well as Potbhare et al , , 

the 1gg/2 expected to play an important role 

vis-^-vis the structure of the yrast levels in Ni isotopes 

with A > 60, However, the non-availability of the core- 


renormalized effective interactions for the (1fyy2t ‘^Pq/ 2 * 
2pi/2» "^^9/2^ space has prevented us from examining 

here the explicit role of the 


We have also presented in Figure III, 2 the results 
obtained by Mooy and Glaudemans^ in the framework of the 

shell model involving the[(2p2/2» ^^1/2* ^^5/2^'^- 
(lf^/2)"^ (2P3/2, 2Pi/2’ configurations. These 

calculations employed the effective charges ( ©y) = (2,60, 

0,72), (3,30, 0,50) and (1,85, 0i82) for the nuclei 
08,60,62^^^ respectively. 
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I 

f 

Considering first the momentuiD-transf er range 

-1 

0,5 < q < 2.1 fm , we find that the projected HFB and 

the shell model estimates are quite close except around 

the second maximum ( q cv, 1,5 fnT^ ) in the nucleus ^^Ni, 

where the latter are smaller by roughly a factor of two. | 

For larger momentum-transf ers the shell model estimates ; 

are an order of magnitude larger — a feature that I 

I 

emphasizes the desirability of precise form factor data I 

for q > 2.1 fm ^ in the nuclei 

It is interesting to note that the proton effective f 

I 

charges employed in the shell model calculations in the ^ 

nuclei are considerably larger than the A-independent , 

value ( e^ = 1 ,8) employed in the present work. This can 

be qualitatively understood in the following manner. The 

shell model value of the proton effective charge simulates 

partly the effects due to the excitations involving two or 

more 1f7/2 "“Photons , Since the proton-neutron inteioction 

are strongest for orbits with large spatial overlap* such 

. 58 , 60 ^. 

processes are expected to be favoured in the nuclei -• ? 

where the energy- difference between the proton and neutron 

Fermi surfaces is not very large » In heavier nuclei 

(A > 62) with larger neutron-excess* the amplitude of proton 

excitations leading to correlated proton-neutron configurations 

is expected to decrease duo to larger gap between the proton 

and neutron Fermi surfaces, ; 
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Figure III, 3 presents a comparison of the calculated 

2 9 j. 

and the experimental lF(q)r values for the 0 -* 2t transi- 

1 

tions in the nuclei These foim factors have 

earlier been analyzed v/ith phenomenological models the 
modified Tassie model and the Gaussian model^ , The 
best-fit calculations carried out in the framework of the 
latter model have also been shown herei 

As in the case of the Ni isotopes, the present form 
factor calculations for the Zn isotopes also employ effective 
charges such that the isovector effective charge, defined by 
(e - e ), is alvrays le. However, it turns out that a 

^ % V 

reasonable variation of the isoscalar effective charge 

with Cj, = ( 0,75, 0,55, 0,75) for the nuclei 

respectively ■ — is required for optimizing the agreement 

with the available B( E2j O"^-* 2^) values. With these 

-1 

effective charges, the form factors for q < 1 ,2 fm computed 
with the self-consistent wave functions are also in excellent 
agreement with the experiments • We notice, however, syste- 
matic discrepancies between the calculated and observed 
form factors for q > 1 «2 fm The calculated form factors 
for q «'* 1 ,5 fm~^ are smaller by about 40 percent than the 
observed ones, and the second maxima of the theoretical 
form factors display small but noticeable shift towards.. 

higher momentum-transfer. These discrepancies may again be 

, - 40- 

reflecting the necessity of explicit Involvement of Ca core- 
excitations as well as the ISg/g configurations i 




FIG. III. 3 EXPERIMENTAL AND CALCULATED FORM FACTORS FOR THE 

o'" TRANSITIONS IN THE NUCLEI > ' Zn. THE 

BROKEN CURVES REPRESENT THE BEST-Fn CALCULATIONS 
WITH THE FOURIER-BESSEL EXPANSION OF THE TRANSITION 
CHARGE DENSITIES. 
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Present calculations reveal (see the second column of 
Table III.1) that the total Intrinsic quadrupole moments 
for the Zn isotopes are substantially larger than the 
ones obtained for the lighter Ni isotopes. This onset 
of sizable oblate deformation may signal enhanced 1gg/.2 
occupation via a lowering of the (downsloping) Nilsson 
orbitals with k = + 9/2, 

As discussed later, the apparent inadequacy of the 
effective charges to simulate the excitations from the 
core into model space - or out of model space into higher 
orbits - in the context of large momentunwtransf er data 
^manifests itself in the form of discrepancies between the 
calculated and the empirical transition densities in the 
interior reg5 on ( r <3 fm) , 

111,4,2 The Reduced Transition Probabilities, B( E2; o"*"** 2^) 

I. 1 I ..- M . . n ., 

and the Static Quadrupole Moments, 0(2^) 

In Table III,1 and Figure III,4 we have presented 
the results for the static quadrupole moments, 0(2^), as 
well as the reduced transition probabilities, B( E2; 0 ** 2^ ) 
in the nuclei 58i60,62j^j 64,66,68^^^ and 

the B(E2) values have been computed wilii the effective charges 
employed in the calculation of the form factors for various 
isotopes discussed in the preceding section. The calculated 
B(E2j 2|) values are seen to be in excellent agreement 



70 


CO o 
c CD 4^ 4-> 
INI XI -H 


fZ 

D 


• C 

CO •H 
C 

o c 

•H CD 
-P > 


ZJ CO 

: c 

c o 

H M 

4-> 
C D 

CD CD 


X5 fp •H 


o c c 

•H iH Q' O 

12 : fO CD -P 

CN O SZ O 

vO p 

CD CL 

O O > 

-P CO CD 

^ o x: o 

00 nj d 

m <H <0 0) 

CD f— I 

•H s D a 

0) p rP > 

»H O 05 

CD <P > O 

0 X 

c 0 ^ -p 

^ x: + ^ 

Cl> -p c>» <p 

x: >^0 

-P C C3' 


•H TJ x: -H 
CD -p 4 ^ 

^ D 

CQ O T5 Xi 

P-. r-H c »H 

n: ou 03 p 

CL. S 4> 

_ 0 ^ c 

T3 rt O 

0 CO £ o 

tH O C|-< 
rH Cn ♦ 0 

0 PCM x: 

Xi 05 0 4 ^ 

0 x: ^ 

r-4 o (O 

q-f 0 
0 O > 

(f) > -H 

0 -H 0 Cn 
D -P 4J 

f-H O •H -''^N 

0 0 C 

> <p D 
tp 

^ 0 CCM O 


CM <H 

O 0 


H O' 

m V 


P 

4-> 0 


1 - CD ^ 

O 0 (0 ^ 

0CM O 

•^0^00 
CM x: rp x> 

jXl 4 J fp V 05 


CQ x; 
“P 


H -f T- P O 
5 : CM CD »H 


0 

P O 


^ 0 ^ 

Hh T~ 4-^ -f* 
CM 0 O 


^ i2 
c 


CJ D 


0 o ^ x> 

CM^ C 
D " 6 'H 

O (N 4-1 
H W #0 

0 0 x: 


H O CQ 



Reference 8 } Reference 9j neference 18 



(Ulo 01 •)(,Z-P'23)a 

Qq^ 6 + 4- y 



fig.iii.4 graphical presentation of the results for electric 

QUADRUPOLE TRANSITION PROBABILITIES AND STATIC 
QUADRUPOLE MOMENTS IN SOME DOUBLY EVEN Ni AND Zn 
ISOTOPES. THE STRAIGHT LINES JOIN THE POINTS 
CALCULATED WITH THE EFFECTIVE CHARGES EMPLOYED IN 
THE FORM FACTOR CALCULATIONS . 
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with the experiiflents j the only noticeable discrepancy 
occurs in the case of the nucleus ^^Zn where the present 
calculation underest.irnates the lower bound of the observed 
value by about 6,6 percent. 

The usefulness of the available Q(2|) vaHues - and 
these are just three in nuniber - vis-^~vis a test of the 
present microscopic description is severely constrained by 
the large error bars arising due to the uncertainties 
associated with the effects of higher excited states in the 
Coulomb excitation processes. The calculated values are, 
however, qualitatively consistent with the available data? 
in particular, the latter seem to support our prediction 
of oblate intrinsic shapes in the nuclei 


111,4,3 Transition Charge Densities 

In the preceding sections we have seen that the 
introduction of the state- and (nearly) mass-independent 
effective charges for the valence particles is on the 
whole quite adequate in most of the cases for describing 
the observed B( E2) values and the form factors. The effective 
charge model implicitly assumes that the contributions arising 
from the core-polarizations or the extra-modol-space compo- 
nents are proportional to the valence space contributions, 

A more realistic model involves rhe mixing of the zeroth— 
order valence wave function with the 2 -pole, njilm giant 
resonances^ this description is consistent with the use of 
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the separable, inultipcle-multipole interaction in conjun- 
ction with the first-order perturbation theory. In view 
of the expected sensitivity of the transition densities 
towards the details of the radial characteristics of various 
parts of wave functions, we have also considered here this 
prescription for invoking core-polarization by assuming 
the lassie model for the 2^-pole excitations. The total 
transition density is constructed as the sum of two terms — 
the valence-space termA(r) with isoscalar effective charge 
set to zero, and the Tassie core-polarization transition 
density term B( r) (Ref. 19): 

p(r) = A(r) + N B(r) (III.l) 

when B( r) is related to the ground state charge densities 
Pgs(^) by the relation^’ ^ 

B(r) = Cm.2) 

The normalization constant M is obtained by reguiring 
overall consistency with the usual effective charge model; 
we adjust N so that the total charge resulting from an 
integration of P ( r) is just e^ , 

In Figure III. 5 we present the transition charge 
densities for Ni isotopes calculated by combining the projected 
HFB predictions with one or the other of the two models for 
the core-polarization component. We have compared our results 
with the empirical transition densities extracted from the 




y (fm) 

FIG. III. 5 TRANSITION CHARGE DENSITIES FOR THE FIRST 2 * STATE 
jjj 58,60,62^. calculated WITH THE EFFECTIVE CHARGE 
model (DOT-DASHED LINE) AND WITH THE TASSIE MODEL 
(dashed line), the EMPIRICAL CHARGE DENSITIES HAVE 
BEEN SHOWN IN THE FORM OF AN ERROR BAND. THE DOTTED 
CURVE REPRESENTS THE SHELL MOEEL RESULTS. 
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avBilsblo form factor data in tho framswork of Fourior— 

Bessel method by Heisenbergk It is seen that in the 
region 3 fm < r < 8 fm the predictions of the effective 
charge ( dot— dashed line) and the lassie prescriptions 
(dashed line) do not differ significantly. Whereas the 
latter prescription yields transition densities in excellent 
agreement with the empirical ones, the estimates involving 
the A-independent effective charges display only, minor 
discrepancies in In the interior region (0<r<3 fm), 

however, the lassie model for core-polarization provides 
vastly improved description of the transition densities, in 
keen contrast with the effective charge model estimates which 
show large oscillations. The effective charge model emphasi- 
zes the role of the underlying valence orbits. On the other 
hand, the collective contributions inherent in the Tassie 
model reduce considerably the single-particle features of 
the valence part. 


V/e have also displayed in Figure III, 5 the shell 
model results for transition densities obtained by Mooy 
and Glaudemans^ with the effective charge prescription for 
the core— polarization. The results are in qualitative agree- 
ment with the projected HFB estimates involving the effective 
charges • The quantitative differences arise frcxn the enhanced 
(lf7/2 - 1f5/2f and (1f7/2-* 2P3/2)'' transitions allowed 
in the present calculations. As pointed out by Mooy and 
Glaudemans^, these transitions govern to a large extent the 






FIG. III. 6 THE EMPIRICAL TRANSITION DENSITIES (ERROR BAND) 
AND THE THEORETICAL PREDICTIONS INVOLVING THE 
EFFECTIVE CHARGE MODEL (DO: -DASHED LINFS) AS WELL 
AC THE TASSIE MODEL (DOTTE'' LINES). 
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first maximum at r 1 fm and the first minimum at 
r ~ 2,5 fm in the valence part of the transition densities. 

In Figure III, 6 we have presented a comparison of 
the projected HFB results for the transition densities in 
the nuclei * * Zn with the empirical values extracted 

from the form factor data by Neuhausen^. Here again we 
find that the lassie prescription for the core participation 
permits a significantly improved description of the empirical 
transition densities. 

Conclusions 

We have discussed here the calculation of the form 

“i" 

factors associated with the electroexcitation of the 2^ 

levels in some doubly even Ni and Zn isotopes. The B(E2) 

values as well as the transition charge densities for 

the o"^ -* 2^ transition have also been calculated. It 

turns out that the projected HFB ansatz for the wave functions 

in conjunction with the realistic 2p-1f shell effective 

interactions provides a nearly parameteivfree framework for 

interpreting and correlating the available data in a 

microscopic perspective. The calculations discussed here 

also demonstrate that the projected HFB technique is equally 

efficacious for the nuclei in the first- as well as the 

20 

second-half of the 2p-1f shell • 
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Present calculations have revealed significant 
discrepancies vis~^~vis the large momentunv-transf er 
(q > 2.0 fm ) form factor data — and the related 
transition densities data for the interior region ( r < 3) 
“"■ in some cases. An explicit inclusion of the "^^Ca core- 

V 

excitations as well as (lgg/2^ configurations may lead 
to an improved description of the form factors and the 
transition densities, apart from reducing considerably 
the quantitative significance of the choice of a model 
for "effective charge" contributions. 
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CHAPTER IV 


i^Np_ tra n sition charge densities assoc iated with 

THE .ELECTROEXC ITATIQM OF THE YRAST LEVELS WITH t'R' = a"** tm 
SOME 2p~1f SHELL NUasi — : 

IV, 1 Introduction 

In the context of a general discussion of the 
shape-collective aspects of nuclear dynamics, a study of 
the electric hexadecupole matrix, elements is a logical 
extension of the studies involving the electric guadrupole 
operator. The non-availability of the B(E4; 0*^- 4'*') 
values has, however, hindered progress on this topic in 
light— and medium— mass nuclei^ the E4 decay branches are 
unmeasurably small with respect to the competing E2 
branches , 

1 —3 

A number of recent experiments have yielded 

important information involving electroexcitation form 

factors covering a broad momentum-transfer range, with 
—1 -f + 

q 5= 0, 5-3.0 fm , for the 0 -* 4 transitions. Apart from 

providing an estimate of the heretofore unmeasured E4 

transition probabilities via an extension of the data into 

the momentum-transf er of the photon point, q = E - E , , 

4 0 

the measured form factors are also expected to be extremely 
useful in the context of the test of various microscopic 
models of the yrast J® = 4**" state. 
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In this chapter we discuss the calculation of the 

IFl^ ( O"^ -♦ 4 '*') values in the nuclei 

Fe and Ni with a view to tost the efficacy of the 

proj Gcted HFB description vis-a-vis the recently-measured 

data. The recent v^/ork of Mooy and Glaudemans"^ has revealed 

that a restricted shell model description of the ( -* 4’^) 

52 54 

form factors in the nuclei Cr, Fe and Ni requires 
highly mass-dependent effective charges which are also 
drastically different from the set required for a reasonable 
quantitative agreement vdth the observed (0 2 ) form 

factors in these nuclei; the values [e.j^(o'^-» 2 "^), © 71 ( 0 "^ ■* 4"^) ] 
required to optimize the fit to the observed data are [1.33, 
2.06], [ 1 .42, 1 ,76] and [3.33, ,1.54] for the isotopes '^^Cr, 
^"^Fe and respectively. In this context we find that 

the HFB description for the = 4 states in the Ti, Cr, 

Fe and Mi isotopes not only obviates the necessity of 
invoking significant mass-dependence, it also permits a 
good quantitative discussion of the (0 -♦ 4 ) form factors 
in terms of affective charges quite close to the ones required 
for a comparable fit to the (0 •* 2 ) data in these nuclei. 

The calculations discussed here have also been carried 

out employing the MVJH interaction described in the preceding 

Chapter, 

In Section IV, 2 we discuss some qualitative features 
associated with the O"*" •*■4 transitions in the 2p-1f shell 
nuclei. In Section IV*3 we discuss the 0 4 transition in 
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the 2p-1f shell nuclei with N or Z = 28. The hexadecupole 
electroexcitation of the 2p-lf shell nuclei with open 
proton and neutron shells is next discussed in the Section 
IV. 4, In view of the recent revival of interest in the 
B(E4j 0 4 ) values due to the availability of the 

2 + -i- 

I FI (0 -* 4 ) values, one is prompted to ask if the usual 
Bohr and Mottelson ( BM) prescription for relating the 
B(E2jO 1 ) values and the intrinsic guadrupole moments 
could also be extended to the hexadecupole operator. In 
Section IV. 5 v^e demonstrate the quantitative as. well as 
qualitative unreliability of the BM prescription vis-^-vis 
the hexadecupole operator by carrying out Hartree-Fock- 
Bogoliubov calculations follovA/ed by explicit angular 
momentum projection for a number of 2p-1f shell nuclei. 

Section IV. 6 contains some concluding remarks. 

IV. 2 Qualitatj-ve Featuros Associated with the O'*' 4'*' 
Tradition's in the' '^p~1f Shelf Mucloi 

"f" "f* 

The qualitative aspects of the 0 ■* 4 form factors 

in the 2p~1f shell nuclei are mainly governed by the two 

radial matrix elements < 1fl j 4 ( qr)I1f > and < 1fl j^C qr)f 2p > , 

if we assume that the q-dependence of these matrix elements 

40 

is not significantly changed by the Ca core-polarization 
processes. It turns out that (see Fig. IV. 1) ”tiereas the 
matrix element < 1fl j 4 (qr)l 1f> possesses just one maximum 
at q 1 .2 fm""^, the matrix element < 1fl j 4 ( qr) I2p > possesses 
a maximum at q r. 1 .1 fm"'"’ followed by a minimum at q 2.2 fm~^ . 
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The zeros of the matrix elements < 1f ! qr) I2p> and 

f ^ ^ 1 f ^ occur at q 1 ,8 and 2,4 fm respectively. 

Since the description of F( q) in the projected HFB framework 
involves e linesr cornbination of the matrix elements 

<1f [j4 (qr)l 1 f> and < 1f t j^C qr) l2p> , it is expected to 
display the follovang general features: 

( i) In view of the proximity of the maxima of the two 
matrix elements, we do not expect significant 
A-dependence in the. position of the first maximum^ 

This feature is strikingly reflected in the observed 

1 ~3 

data ; the first maximum appears at practically 

-1 

the same position — at q ~ 1 ,2 fm — in various 
nuclei, 

(ii) We do not anticipate the occurrence of the second 
maximum in nuclei with structure dominated by the 
(if?/^)*^ configurations for protons as well as 
neutrons, since, in such cases, the form factor 

is governed by only the matrix element < 1f Ij4 (qr)l 1 f> , 

(iii) In the deformed 2p-1f shell nuclei, the 2p2^2 
orbit is also expected to play an important role. 

The magnitude of lF(q)I^ at the second maximum as 
well as the location of this maximum is, therefore, 
expected to display considerable variation over the 
range q = 1 .8 - 2.5 fm''^ depending on the relative 
importance of the matrix element < 1fl j4( gr)I 2p > i 
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ligMde^cuioole Electroexcitati nn in the 2n-if shell 
Mucloi with N or Y = — - 

IV, 3,1 Form Factors 

The squared form factors for the 0^ 4'*' transitions 

in the nuclei Ti, ' Cr, Fe and ^^Ni are calculated and 
compared with the experiments^”^ in Fig, IV, 2, In the 
calculations reported in, this and the next section, we have 
employed the mass-independent set of effective charges 
(e^, e^) = (1,5, 0,5) for the nuclei and 50,52,54^^^ 

In the case of the nuclei ^"^Fe and ^*^Ni, the values employed 
are (1,7, 0,7) and (1,8, 0,8), respectively. 

From the results presented in Figure IV, 2 one finds 

that the calculated form factors are in very good quantita- 

2 

tivG agreement with the experiments . In the case of the 
52 

nucleus Cr present calculation predicts correctly the 
position as well as the magnitude of the second maximum. 

It is seen that a lack of enough reliable data points 
for q < 1 ,0 fm"^ prevents a detailed test of the calculated 

“f" “f- 

results. The importance of the lov\^-q data for the 0 4 

transitions is further emphasized by the fact that the 
experimental values for B(E4; 0 •* 4 ) and these are the 
limits of the relevant form factors as q approaches zero -*- 
are not available for testing the microscopic description. 
This contrasts keenly with the situation in the case of the 
o'" -* 2*" transitions; in these cases the lF(q)I^ estimates 
for q < 1 are usually supplemented by the B(E2) values 
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We find that in the nuclei considered here the 
magnitude cf the calculated IF(q)I^ at the second maximum 
is smaller by more than two orders of magnitude compared 
to the values at the first maximum. In view of the general 
discussion of the form factors given earlier, the result 
implies the dominance of the configurations in 

the yrast states with = 0 , 4"^ in these nuclei. 

We have also given in Figure IV, 2 the results 

obtained by Mooy and Glaudemans"^ in the cases of the nuclei 
52 54 60 

Cr, Fe and Ni, The effective charges ( = (2,06, 

(1,76, 0,76) and(l,54, 0^54) were employed for 
52 54 60 >. 

the nuclei Cr, Fe and Ni, respectively, in these 
calculations. The shell model calculations have been 
reported only for the range q = 0,6 - 1 ,8 fm covering 
just the first maxima. The calculations reproduce the 
observed data quite well in this momentum-transf er range. 

We note that in the nucleus Cr, the shell model (e^, 
values are considerably larger than the ones considered in 
the present calculation. Largo renormalized values of the 
effective charge, which imply an enhanced involvement of 
the 2 p 2^2 orbit, are not unexpected in this nucleus in view 
of the significant non-sphericity of the zeroth-order 
description of the proton wave function in terms of the 

/ 41^ 

( 1 f7/2) configuration. 
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It may be pointed out here that the effective 
charges employed in the present calculation of the 4''’ 
form factors for various 2p~lf shell nuclei are either 
the same or lie within 12 percent of the values employed 
in our earlier calculations'^ for the 0 ^ •* 2 ^ form factors 
in these nuclei. This feature is not shared by the shell 
model calculations (see columns 5 and 7, Table IV, l); 
tho shell model results are characterized by large differ- 
ences between the effective charges employed for the 
** 1 “ 

0 -* 2 transitions and those for the 0 4 transition, 

IV, 3, 2 Transition Charge Densities 

In Figure IV,3 we have displayed the calculated 

•> 1 “ 

transition charge densities associated with the 0 -* 4 

50 52 54 60 

transitions in the nuclei ' Ti, Cr,. Fe and Mi, The 
transition charge densities have been computed using the 
effective charges employed in the form factor calculations 

In view of the nonavailability of tho empirical 

0"^^. 4 ^' transition densities extracted from the Fourier- 

Bessel analyses, it is not possible to assess the results 

obtained here. However, it is interesting to note that the 

present calculation does reveal significant qualitative 

differences betv^een the results obtained for the tradi- 

50 52 54-, 

tional nuclei — the nuclei Ti, Cr and Te — 

and those for tho nucleus ^°Ni; the presence of tho 
configurations in the latter nucleus is signalled by the 
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appearance of an oscillatory form. 


IV, 4 Hexadecupole Electroexcitation in the Nuclei 


^°Cr and^Cr 


IV, 4,1 Form Factors 


Figure IV, 4 compares the calculated and observed IF I' 

"f" "f" “^8 

for the 0 -* 4 transitions in the nuclei Ti, Cr and 


-f + 4R 

Although the 0 •* 2 transitions in the nuclei Ti 
and ^^Cr have earlier been discussed by Iwaraoto et al ,*^ 
in the framework of the shell model involying restricted 
configurations, no shell model calculation has so far been 

"4* “f* 

reported in these nuclei for the 0 -* 4 form factors. 

Since these nuclei do not involve the 1 ^ 7/2 sub~shell 
closure, it is difficult to justify a priori any truncation 
scheme for reducing the number of shell model configurationsi 


As mentioned earlier, the projected HFB calculations 
have been carried out with a constant set of effective 
charges with e^ =1.5 and = 0,5 — for the nuclei 

considered here. It is seen that the present calculation 
is quite successful in reproducing the magnitude of the 
form factors at the first maxima. However, one observes 

-1 

discrepancies in the momentum-transfer range 1 ,5 < q <2,7 fm ; 
the position of the first minimum is shifted towards hi^er 
momentum-transfers and the calculation underestimates the 
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2 

magnitude of Ipl around its second maximum in the nuclei 
5 0 54 

’ Cr, It is seen that the quadrupole momaits of the 
intrinsic states (see column 2, Table IV, 1) in the nuclei 
"^^i, are significantly larger than those for the 

nuclei with N or Z = 28, The onset of sizable quadiupole 
deformations may warrant the inclusion of the 199/2 
in the configuration space. It would also be worthwhile 
to examine the effect of incorporating more configurational 
admixtures within 2p-1f shell by considering quasiparticle- 
excited K =! 0 intrinsic states, in addition to the usual 
ones resulting from the HFB calculation, 

IV ,4. 2 Transition Charge Densities 

In Figure IV,5 we have displayed the calculated 

+ + 

transition charge densities associated with the 0 -* 4 
transition in the nuclei *^^Ti and 

One notices in these nuclei significant correlation 

between the ratio of the form factors it their first and 

second maxima and the size of the ’dip’ in the calculated 

transition densities in the interior region (0 <r< 3 fm); 

the near-disappearance of the dip (at r^ 2 fm) in the 

50 

transition charge density in Cr is seen to be associated 
with the significantly reduced estimates (relative to the 
values around first maximum) of the form factors around 
second maximum. The empirical transition charge densities 
are thus expected to provide important clues regarding the 



J>(efm^) X 10 


1.0 
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origin of the discrepancy between the calculated and 

-1 

the observed form factors for q > 1 ,5 fm in the nucleus 



The rotational spectra in nuclei are usually 

discussed in the framework of the unified model of Bohr 
1 0 

and Mottelson , The model is characterized by the 
decoupling of single- particle and rotational motions; 
the wave function of a yrast state in quasirotational 
nuclei is given, for example, by the product of the 
symmetric top v\ave function and an intrinsic state descri- 
bing, respectively, the collective rotation and the single- 
particle motion. The BM or the rotor model relates the 
reduced electric 2^-pole transition probabilities to the 
intrinsic 2^-pole moments ““ the expectation values of 
operators with respect to the intrinsic states ■— in a 
simple manners 

^ _2 
Ji L 

B(EL; jt jJ) = (1/16-rt) I 

Lo 0 0 J 

(IV.1) 

1 1 

As discussed by Vi liars , eq. (IV.I) can also be 
realized in a completely microscopic (quantal) perspective 
by considering^ the zeroth-order approximation to the matrix 
element < § 1 P ^ 1^ > in powers of < ^ IJ^l § > , 

where ] S > is an intrinsic state and P'^ is the Peierls-Yoccoz 
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projection operator. It nay be stressed here that the 

apparent validity of this approximation scheme depends 
2 

obly on the <J > value and not on the multipole involved. 

The relation ( IV.l ) has been extensively used in the 

1 9 

past for obtaining semiquantitative estimates of the 

reduced transition probabilities for electric quadrupole 

transitions, B(E 25 o”^ 2'"), in terms of intrinsic quad~ 

2 

rupole moments, <Q^>, resulting fron either phenomenological 
models such as the Nilsson model or microscopic, Hartree- 
Fock (HF) or Hartree-Fock-Bogoliubov (HFB) descriptions. 

The BM prescription has also often been invoked for 
obtaining Nilsson parameters for quadrupole deformations 
in' different mass-regions by first extracting the <5 IQ^ I§> 
values from the measured transition probabilities , 

As mentioned in the beginning of this Chapter, the 
recent accumulation of the 0 •* 4 form factor data for 
a number of 2p-1f shell nuclei implies, in a somGv\tiat 
indirect manner, the availability of the heretofore 
unmeasured E4 transition probabilities, B(E4; o"*"-* 4"^), 
since the latter can be looked upon as the appropriate 
q 0 limits of the IF! Co' -* 4 ) values. In this context 
one is tempted to enquire if the equation (IV,t) could also 
be used for relating the E4 reduced transition probabilities 
and the hexadecupole moments of the intrinsic states. 

In this section we present explicit calculations to 
show that, viiereas the zeroth-order approximation implicit 
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in the derivation of eq, (IV«1) in the Peierls-Yoccoz 

1 3 

projection formalism is qualitatively as well as, to 

a good extent, quantitatively quite reliable in the case.. 

of the quadrupole operator, it 'breaks down completely in 

the case of the hexadecupole operator, irrespective of 
2 

the value of < J > , It is pointed out that the validity 
of various orders of approximation depends primarily on 
the s imilarities between a pure quadrupole field and 
the (one-body) self-consistent HF field resulting from 
the effective interaction and not so much on the magnitude 
of the < value. The anomalous nature of the results for 
the hexadecupole operator in the 2p>-1 f shell can be traced 
back to the significant di ff erenc es between the energy 
sequence as well as the structure of the eigenstates of 
the Qq and fields in the shell. 


Vife have listed in Table IV, 2 the intrinsic quadrupole 
as well as hexadecupo.le moments of the axially symmetric 
HFB states (resulting from the Will interaction) for the 
nuclei ’ ’ Ti, ’ ’ Cr and * Fe. We have also 

presented the reduced E2 and E4 transition probabilities 
resulting from explici-t angular momentum projection, as 
well as their estimates resulting from the BM prescription 
( eq. ( IV,1 ) ) , 


A graphical presentation of the results (see Figure .. 
IV, 6) shov-ys that the rotor model predictions are in excellent 
qualitative agreement with the PI-iFB results fcBC the 



TABLE IV. 2 


1 00 


The intrinsic multipole moments as well as the reduced 

transition probabilities involving the quadrupole and the 

hexadecupole operators for some doubly even 2p-1f shell nuclei. 

The • reduced transition probabilities resulting from explicit 

angular momentum projection on the HFB intrinsic states have 

been compared with the rotor model predictions. The intrinsic 

2^-pole moments have been given in units of b^, where 

b = iYi/mcS) ' is the oscillator parameter. The reduced transition 

probabilities B( E2| O"*" -*■ 2**) and B(E4j O"^ -*■ 4"^) have been given 
O —‘^0 A 9 —109 R 

in units of e^.10 cm and e ,10 cm , respectively 


Nuc leus 

9 

<<> 

B( E2; 

Projected 

HFB 

o'*"- 2 *) 

Rotor 

Model 

A 

JO 

O 4^ 
V 

B(E4;0'^-* 

Projected 

HFB 

4-^) 

Rotor 

Model 


16.4 

4,7 

3.7 

21 .3 

9.6 

156.0 

5°T1 

-9,4 

2.5 

1 .2 

18.9 

25.6 

79.2 


15,7 

4.9 

3,6 

20.0 

14.9 

150i4 

50, 

Cr 

26.3 

10.7 

9.7 

23.2 

5.2 

192.2 

52 

Cr 

11.4 

3.7 

1 .9 

15 ,7 

13,7 

92.6 


26.7 

12.5 

10.6 

33.7 

13.8 

449.2 


20.8 

8.0 

6.2 

-19.3 

13.1 

140.5 


9.0 

3.0 

1 .2 

-7.5 

14.1 

22.3 




B(E4 i ) («2 cm®) B(E 2 jOW) (n^ iS^^cm*^) 



FIG. IV. 6 COMPARISON OF THE REDUCED TRANSITION PROBABILITIES 
RESULTING FROM THE PHFB AS WELL AS THE ROTOR MODELS. 
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B( E 25 0 “* 2 ) values^ In fact, for the N ^ 28 nuclei, 
even the quantitative agreement between the rotor model 
and the PHFB estimates is quite goodj the maximum discre- 
pancy is only about 27 percent of the latter. Since the 
nuclei with N 28 also possess significantly larger 
intrinsic quadrupole moment — and consequently larger <J% 
values compared to the nuclei with = 28, the results., 
obtained here are consistent v/ith our h priori expectation 
concerning the general validity of the rotor model. 

In sharp contrast to the situation vis-a-vis the 

B( E2j 0 *• 2 ) values, we observe (see Figure IV, 6) dramatic 

failure both in the quantitative as well as qualitative 

4 

sense — of the BM prescription as applied to the 
operator. In fact the rotor model estimates are larger by 
an order of magnitude than the PHFB estimates in most of 
the cases. 

In v'\hat follo’vs v;g make an attempt to examine 
qualitatively the reasons for the inadequacy of the rotor 
model prescription as applied to the hexadecupole operator. 
In the framework of the rotor model, the state of angular 
momentum J and projection U (in the laboratory frame) is 
given by 

l*MO>= [(2J+1)/16n2] CIV,2) 

Using these wave functions the matrix elements of theQ^ 

*}T/ -j- jfD -|- 

operator between = 0 and L state becomes 
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< 


L 

00 



1^0 > 
00 


"l L o" 
o' 0 0 


< 5 I qL ] § > 

0 0 o 


( IV.3) 


We next employ the angular momentum projected states 

i A 

given by equation (11,128). As shown by Ripka , the matrix 
elements of the operator between projected states are 
given by 


<- iQ^ ! > 

^ ^00 'o' 00 ^ 




•1 /2 ^ 

/ dp 

1^ o 

-ipj. 


T — T 

Slnp {P)<sy e y chis > 


li- 0 
(IV.4} 


\Ahere 


ar = 


r d ' ^ T -iP J . , 

dp SinP d^ (P ) <^^Ie YI ^ > 

2. ^ 00 ^ ' o o 


(IV.5) 


A comparison of equaxions ( IV,C) and ( IV, 4) shows 
that they are expected to yield nearly the same results in 
the zeroth-order approximation and this involves the usual 
assumption that the integrals give nonvanishing contribution 
only for p»^0 provided the intrinsic states 1 ^ q> are 

eigenstates ( or near-eigenstates) of the operator , 

(Note that in the zeroth-order approximation the p -summation 
is also automatically restricted to just the p= 0 term). 

In Table IV, 3 we have displayed the structure of the single^; 
particle states in the cuadrupole and hexadecupole fields in 

„4 

the 2p-1f shell. It is seen that the eigenstates of the 
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TABLE IV, 3 

Structure of single- particle states in pure quadrupole and 
hexadecupole fields in the 2p-1f shell,- Here k=<l^> is 
the expectation value of the component of the angular 
momentum along the symmetry axis and ^ 2 ^ ^4^ eigen- 

value of the intrinsic quadrupole (hexadecupole) moments in 
units of b^( 


Quadrupole field Hexadecupole field 


k 

"2 

^ ( k, ^2) 

Ic 

®4 

§(k,e4) . 

0 

6*0 

\ 

0.63fQ-0.77PQ 

0 

19.9 

0.80f^-0,59pQ 

±1 

3,0.. 

0o89f^^-0,45p.,^^ 

±1 

9.8 

0.74f^ -0.68p^^ 

mum IbMM 

+2 

0 

• 

0 


+3 

4,5 

00 

0 

0 

• 

0 

0.77fQ+0.63pQ 

+1 

-8.3 

0.68f^^+ p.74p_^^ 

+1 

-3.0 ^ 

0.45f_^,^+0,89p_^^ 

+2 

-10*5 

^+2 

+3 

-3.0 


0 

-10*9 

0,59fQ+ 0.80p^ 
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2 

operator differ substantially from those of the 

operator. Thus the intrinsic states I which are near- 

2 

eigenstates of the operator due to the well-known 

15 

quadrupole-quadrupole dominance of the effective nuclear 
interactions, are not in general the eigenstates of the 

4 

operator. In other words, whereas the nature of the 

interaction facilitates the similarity between the rotor 

model and PHFB predictions for the quadrupole operator (by 

yielding intrinsic states that are close to the eigen- 
2 

states of Qq) , the inherent differences in the 2p-1f shell 

2 4 

eigenstates of and operator do not permit the recovery 
of the rotor limit from the expression (IV.4), 

IV, 6 Conclusions 

We have discussed in this Chapter the calculation of 
the form factors as v.;ell as the transition charge densities 
associated with the olectrooxcitation of the yrast 4"*" 
levels in some doubly even 2p-rf shell nuclei. It turns 
out that the use cf the j’’^ = o"*", 2*, 4‘ states projected 

from variational intrinsic states resulting from realistic 
2p-1f shell intera.ctions permits a reasonably satisfactory 
simultaneous description of both the 0‘ >* 2 as well as 

I ^ 

0 4 form factor data with nearly the same set of effective 


charges 
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The present work has eraphasi 2 ed the necessity of 

providing an improved description of the large-Q part 

+ + 50 54 

of the 0 •* 4 form factors in the nuclei * Cr, In 

this context, the empirical transition charge densities are 

expected to provide valuable guidelines. 
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CHAPTER V 


CONCLUSIONS 


In recent y^rs the experimental studies involving 
electroexcitation of nuclear levels have provided extremely 
valuable information about several heretofore unexamined 
aspects of the nuclear structure of medium-mass nuclei. 

In contrast to the conventional experimental studies 
involving Y-transitions, the momentum- transfer dependence 
of the matrix element ( extracted fran the measured cross 
sections) yields detailed information vis-^-vis the spatial 
structure of tho nuclear states participating in the 
transition. 

Several earlier theoretical investigations have 
involved a comparison of the observed data v/ith the predi- 
ctions of nuclear collective models. The collective model 
analyses have sought to reproduce eloctroexcitation form 
factors by postulating general functional forms of the 
nuclear charge Qtnsities, While condensing the experimental 
data into concise parametrizations , such analyses provide 
us with only a qualitative characterization of the nuclear 
excitations; they do not permit us to interpret and correlate 
the data from the point of view of a '.-.ider perspective 
involving other microscopic features of nuclear stincture. 

In this thesis we have attempted an assessment of 
the capability of the projected HFB technique vis-S-vis the 



109 


interpretation of such fundamental results as those embodied 

in the data on el'-^ctroexcitatipn transition charge densities 

as- well as form factors involving several 2p-1f shell nuclei. 

It turns out that the available form factor data out to about 
-1 

2,5 fm can be reproduced in most of. the cases in a fairl-y 
satisfactory nanner in terms of reasonable values of nearly 
mass -independent effective charges. The calculations presented 
in this thesis have also brought out the decisive role that 
the empirical transition charge densities in the Ni and Zn 
isotopes, extracted from the cfeta via the Fourier-Bessel 
analysis, play vis-h-vis the choice of a model of core- 
polarization contributions. 

The results presented in this thesis have revealed 

significant discrepancies in some cases between the observed 

and the calculated form factors for the large-q regions. 

These discrepancies seem to ’.varrant a careful study of the 

40 

role played by an explicit inclusion of the Ca core- 
excitations as 'Well as th ? configurations involving the l9g/2 
orb-it. Further, it ap-'m^ars u'oethwhilo to attempt a systematic 
renormalization of the transition charge density operator in 
terms of various core-excitation graphs. Apart from obviating 
the dependence on the existing 'models' for the coro-polarizatior; 
transition density, such a study is expected to reveal the 
relative importance of various core-excitation processes ! 

vis-S-vis the detailed coordinate-space structure of the 
transition charge densities. 



APPENDIX A 


THE ^P3/2> ^^1/2* MTRIX ELH AENTS 

OF THE KB VIELL AS THE MH INI ENACTIONS 


The matrix elements <abJTi VI cdJT> of the Kuo- 
Brown (KB) interaction are tabulated* The columns labelled 
by (a,b,c,d) contain twice the j-value of the various 
shell model orbits. 


T a 
1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 

1 7 


b 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

5 

5 


c d 
7 7 

7 5 

7 3 

7 1 

5 5 

5 3 

5 1 

3 3 

3 1 

1 1 
7 5 

7 3 

1 


J 

<IV i > 

J 

<IVl> 

J 

0 

-1 .807 

2 

-0,785 

4 

6 

0.226 




2 

0.000 

4 

-0.406 

6 

2 

-0,502 

4 

-0.307 


4 

-0.293 




0 

-2,788 

2 

-0,638 

4 

2 

0.305 

4 

0,305 


2 

-0i470 




0 

-0.783 

2 

-0.269 


2 

-0.250 




0 

• 

o 

1 




1 

-0,287 

2 

-0.096 

3 

4 

0,031 

. 5 

0.156 

6 

2 

0.014 

3 

-0,118 

4 

5 

-OiOlO 




3 

-0i012 

4 

-0.093 



<l Vl> 
.0i087 

-0*716 

-0i400 


0;022 

-0*894 

-0.114 


1 


7 


5 


7 
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T 

a 

b 

c 

d 

J <IVl> 

J 

<1V1> 

J 

<ivi> 

1 

7 

5 

5 

5 

2 -0;531 

4 

-0.468 



1 

7 

5 

5 

3 

1 -0.035 

4 0.493 

0 

0.295 

3 

-0.132 

1 

7 

5 

5 

1 

2 -0,342 

3 

0.080 



1 

7 

5 

3 

3 

2 -0i018 





1 

7 

5 

3 

1 

1 -0.079 

2 

-0.150 



1 

7 

3 

7 

3 

2 -0.861 

5 0i145 

3 

-0.027 

4 

-0.047 

1 

7 

3 

7 

1 

3 0,078 

4 

-0.502 



1 

7 

3 

5 

5 

2 -0.446 

4 

-0.156 



1 

7 

3 

5 

3 

2 0.402 

3 

0,026 

4 

0.619 

1 

7 

3 

5 

1 

2 -0i966 

3 

-0.059 



1 

7 

3 

3 

3 

2 -0i325 





1 

7 

3 

3 

1 

2 -0.319 





1 

7 

1 

7 

1 

3 0.029 

4 

-0.274 



1 

7 

1 

5 

5 

4 -0.252 





1 

7 

1 

5 

3 

3 0.140 

4 

0.531 



1 

7 

1 

5 

1 

3 -0,078 





1 

5 

5 

5 

5 

0 -O486O 

2 

-0.218 

4 

0.298 

1 

5 

5 

5 

3 

2 0i023 

4 

0.102 



1 

5 

5 

5 

1 

2 -Oil 86 





1 

5 

5 

3 

3 

0 -0i777 

2 

-0.128 



1 

5 

5 

3 

1 

2 -0,256 





1 

5 

5 

1 

1 

0 -0.392 
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T 

a 

b 

c 

d 

J < iv 1 > 

J 

< 1 V 1 > 

J 

<IVI> 

1 

5 

3 

5 

3 

1 -0.032 

4 -0.243 

2 

0,215 

3 

0,120 

1 

5 

3 

5 

1 

2 0i389 

3 

0.007 



1 

5 

3 

3 

3 

2 0i068 





1 

5 

3 

3 

1 

1 -0i058 

2 

0.159 



1 

5 

1 

5 

1 

2 -0,135 

3 

0i205 



1 

5 

1 

3 

3 

2 -0.047 





1 

5 

1 

3 

1 

2 -0.241 





1 

3 

3 

3 

3 

0 -1 .206 

2 

-0,380 



1 

3 

3 

3 

1 

2 -0i601 





1 

3 

3 

1 

1 

0 -1 i465 





1 

3 

1 

3 

1 

1 0.152 

2 

-0,688 



1 

1 

1 

1 

1 

0 -0.249 





0 

7 

7 

7 

7 

1 -0.525 

7 -2.199 

3 

-0.208 

5 

-0.502 

0 

7 

7 

7 

5 

1 1 .894 

3 

1 .005 

5 

0.901 

0 

7 

7 

7 

3 

3 -0i482 

5 

-0.816 



0 

7 

7 

7 

1 

3 0.640 





0 

7 

7 

5 

5 

1 1 .071 

3 

0,517 

5 

0;170 

0 

7 

7 

p; 

u 

3 

1 0.227 

3 

-0.087 



0 

7 

7 

5 

1 

3 0.033 





0 

7 

7 

3 

3 

1 -0i279 

3 

-0.296 



0 

7 

7 

3 

1 

1 0,392 





0 

7 

7 

1 

1 

1 0.184 
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T 

a 

b 

c 

d 

J 

< IVI > 

J 

< IVI > 

J 

< 1 V 1 > 

0 

7 

5 

7 

5 

1 

-3,621 

2 

-2,731 

3 

-0.985 






4 

-1 .886 

5 

-0,112 

6 

-2.217 

0 

7 

5 

7 

3 

2 

-0,735 

3 

0.328 

4 

-0.039 






5 

0.405 





0 

7 

5 

7 

1 

3 

-0,440 

4 

-0.674 



0 

7 

5 

5 

5 

1 

-0.416 

3 

0.675 

5 

1 .138 

0 

7 

5 

5 

3 

1 

0,953 

2 

-0.705 

3 

0.528 






4 

-0.544 





0 

7 

5 

5 

1 

2 

0,427 

3 

0.539 



0 

7 

5 

3 

3 

1 

0,871 

3 

0.538 



0 

7 ■ 

5 

3 

1 

1 

-1 .449 

2 

-0.659 



0 

7 

5 

1 

1 

1 

0,168 





0 

7 

3 

7 

3 

2 

-0,293 

3 

-0.604 

4 

-0.164 






5 

-2.165 





0 

7 

3 

7 

1 

3 

1 .295 

4 

0.108 



0 

7 

3 

5 

5 

3 

0.158 

5 

0.048 



0 

7 

3 

5 

3 

2 

-1 .123 

3 

-0.408 

4 

-0i614 

0 

7 

3 

5 

1 

2 

0.838 

3 

0.518 



0 

7 

3 

3 

3 

3 

-0.511 





0 

7 

3 

3 

1 

2 

-0i481 





0 

7 

1 

7 

1 

3 

-1 .484 

4 

-0.746 



0 

7 

1 

5 

5 

3 

-0.270 





0 

7 

1 

5 

3 

3 

0.139 

4 

-1 .188 



0 

7 

1 

5 

1 

3 

0.400 





0 

7 

1 

3 

3 

3 

0,461 
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T 

a 

b 

c 

d 

J 

< IVI> 

J 

<IVI > 

J 

<ivi> 

0 

5 

5 

5 

5 

1 

-0i012 

3 

-0.196 

5 

-1 i687 

0 

5 

5 

5 

3 

1 

-0»292 

3 

-0,327 



0 

5 

5 

5 

1 

3 

-0.391 





0 

5 

5 

3 

3 

1 

0i040 

3 

-0.182 



0 

5 

5 

3 

1 

1 

0i159 





0 

5 

5 

1 

1 

1 

-0;084 





0 

5 

3 

5 

3 

1 

-2.164 

2 

-1 .236 

3 

-0i459 






4 

-0,722 





0 

5 

3 

5 

1 

2 

0.255 

3 

-0,988 



0 

5 

3 

3 

3 

1 

-0i073 

3 

-0.340 



0 

5 

3 

3 

1 

1 

0.620 

2 

-0,071 



0 

5 

3 

1 

1 

1 

-0i393 





0 

5 

1 

5 

1 

2 

-0.081 

3 

-1.234 



0 

5 

1 

3 

3 

3 

-0.084 





0 

5 

1 

3 

1 

2 

0.464 





0 

3 

3 

3 

3 

1 

-0,639 

3 

-1 .832 



0 

3 

3 

3 

1 

1 

1 i554 





0 

3 

3 

1 

1 

1 

0,709 





0 

3 

1 

3 

1 

1 

-2i291 

2 

-1 .897 



0 

3 

1 

1 

1 

1 

0;686 





0 

1 

1 

1 

1 

1 

-1 .074 
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The MVJH interaction referred to in the thesis is 
obtained from the KB interaction by substituting the 
following matrix elements in place of their counterparts: 


T 

a 

b 

c 

d 

J 

< I V 

J 

<IV1 > 

J 

<ivi> 

1 

7 

7 

7 

7 

0 

-2.107 

2 

-1 ,115 

4 

-0i097 

1 

7 

5 

7 

5 

1 

-0.037 

2 

0,154 

3 

0i272 





4 

0.281 

5 

0.406 

6 

-0.644 

1 

7 

3 

7 

3 

2 

-0;561 

3 

0,253 

4 

0.283 





5 

0i485 





1 

7 

1 

7 

1 

3 

0,279 

4 

-0,024 
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